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within  the  WES  reservation  at  Vicksburg,  Mississippi.  A  number  of  in  situ  direct 
shear  tests  were  conducted  at  each  test  site  and  the  results  used  to  fit  a 
mathematical  soil  model  with  a  rate-dependent  nonlinear  failure  envelope  for 
input  to  the  tracked  vehicle  prediction  model.  Thirty-five  circular-turn  tests 
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from  these  16  "good  tests"  were  used  to  validate  both  the  steady-state  and  the 
transient  versions  of  the  model. 
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For  the  transient  predictions',  Fourier  series  were  fit  to  the  entire 
measured  time  histories  of  the  inner  and  outer  track  velocities  and  used  as  input 
to  the  model.  Calculated  time  histories  of  vehicle  speed,  lateral  acceleration, 
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maximum  vehicle  velocity  based  on  either  specified  stability  criteria  or  power 
available  at  the  sprockets.  Predictions  for  both  transient  and  steady-state 
conditions  compared  very  favorably  with  the  test  data,  indicating  that  the  model 
is  capable  of  predicting  the  steering  performance  of  tracked  vehicles  on  soft 
soil. 
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PREFACE 


This  investigation  was  conducted  for  the  Office,  Chief  of 
Engineers,  U.  S.  Army,  by  personnel  of  the  U.  S.  Army  Engineer 
Waterways  Experiment  Station  (WES),  CE,  as  a  part  of  Project 
4A161102AT22,  Task  CO,  Work  Unit  001,  "Dynamic  Soil-Track  Interactions 
Governing  High-Speed  Tracked  Vehicle  Performance." 

The  mathematical  model,  prediction  methodology,  and  analyses 
reported  herein  were  performed  by  Drs.  Behzad  Rohani  and  George  Y. 
Baladi  of  the  Geomechanics  Division  (GD) ,  Structures  Laboratory,  during 
the  period  October  1980  -  October  1981  under  the  general  direction  of 
Mr.  C.  J.  Nuttall,  Jr.,  Chief,  Mobility  Systems  Division  (MSD) , 
Geotechnical  Laboratory  (GL) ,  and  Dr.  W.  F.  Marcuson  III,  Chief,  GL. 

The  field  test  program  was  directed  by  Messrs.  N.  R.  Murphy,  Jr.,  B.  G. 
Schreiner,  and  C.  E.  Green,  MSD.  The  field  direct  shear  device 
described  in  Appendix  B  for  measurements  of  soil  properties  was 
designed  by  Mr.  J.  Q.  Ehrgott,  GD.  The  field  measurements  of  the 
vehicle  performance  were  processed  by  Mr.  P.  J.  Kuykendall,  MSD. 
Numerical  calculations  using  the  WES  terrain-vehicle  interaction  model 
were  performed  by  Mr.  D.  E.  Barnes  and  Mrs.  J.  T.  Carlisle,  GD.  This 
report  was  written  by  Drs.  Rohani  and  Baladi. 

COL  Nelson  P.  Conover,  CE,  and  COL  Tilf  ,  .  C.  C'.-el,  CE,  were 
Commanders  and  Directors  of  the  WES  during  the  investigation. 

Mr.  Fred  R.  Brown  was  Technical  Director. 
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CONVERSION  FACTORS,  U.  S.  CUSTOMARY  TO  METRIC  (SI) 
UNITS  OF  MEASUREMENT 


U.  S.  customary  units  of  measurement  used  in  this  report  can  be  con¬ 
verted  to  metric  (SI)  units  as  follows: 

_ Multiply _ By _ _ To  Obtain _ 


feet 

horsepower 

inches 

inches  per  second 

miles  per  hour  (U.  S.  statute) 

pounds  (force) 

pounds  (force) -inch-second 
squared 

pounds  (force)  per  cubic  inch 
pounds  (force)  per  square  inch 
pounds  (mass) 

pounds  (mass)  per  cubic  foot 
square  inches 


0.3048 

745.6999 

25.4 

25.4 

1.609344 

4.448222 

0.11306064 

0.2714 

6.894757 

0.4535924 

16.01846 

6.4516 


metres 

watts 

millimetres 

millimetres  per  second 
kilometres  per  hour 
newtons 

kilograms-square  metres 

megapascals  per  metre 

kilopascals 

kilograms 

kilograms  per  cubic  metre 
square  centimetres 


STEERABILITY  ANALYSIS  OF  TRACKED  VEHICLES  ON  SOFT  SOIL; 


THEORETICAL  PREDICTIONS  VERSUS  FIELD  MEASUREMENTS 


PART  I:  INTRODUCTION 


Background 

1.  Development  of  high-mobility/agility  tracked  combat  vehicles 
has  received  considerable  attention  recently  because  of  the  possibili¬ 
ties  these  vehicles  offer  for  increased  battlefield  survivability 
through  the  avoidance,  by  high-speed  and  violent  maneuver,  of  hits  by 
high-velocity  projectiles  and  missiles.  In  order  to  design  and  develop 
such  vehicles  rationally,  it  is  necessary  to  have  a  quantitative  under¬ 
standing  of  the  interrelationship  between  the  terrain  factors  (soil 
type,  soil  shear  strength,  and  compressibility,  etc.)  and  the  vehicle 
characteristics  (weight,  track  length  and  width,  location  of  center  of 
gravity,  etc.)  during  steering.  The  actual  mechanism  of  terrain- 
vehicle  interaction  during  steering  is  undoubtedly  very  complex.  Thus, 
in  order  to  study  such  an  interrelationship,  it  is  necessary  to  con¬ 
struct  idealized  mathematical  models  of  the  actual  system.  The  accuracy 
and  range  of  application  of  such  models  must,  of  course,  be  determined 
from  actual  mobility  experiments  and  obviously  must  depend  on  the 
degree  of  relevance  of  the  idealized  model  as  an  approximation  to  the 
real  behavior.  A  research  program  was  initiated  at  the  U.  S.  Army 
Engineer  Waterways  Experiment  Station  (WES)  in  1976  to  develop  a  mathe¬ 
matical  model  of  terrain-vehicle  interaction  for  predicting  the  steering 
performance  of  ground-crawling  vehicles  operating  off  the  road.  Devel¬ 
opment  of  the  model  was  completed  in  1978  (Baladi  and  Rohani,  1979). 

For  its  initial  validation,  results  from  only  five  circular-turn  tests, 
all  conducted  at  one  site,  were  available  for  comparison  with  theoreti¬ 
cal  predictions  (Baladi  and  Rohani,  1981).  Data  from  35  tests  conducted 
at  three  different  soft-soil  sites  have  recently  been  reduced  and  are 
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now  available  for  further  investigation  of  the  model's  validity 
(Schreiner  and  Green,  1980;  Green  (in  preparation)). 


Purpose  and  Scope 


2.  The  purpose  of  this  report  is  to  assess  the  ability  of  the  WES 
terrain-vehicle  interaction  model  to  predict  the  steering  performance 
of  tracked  vehicles  on  soft  soils  by  comparing  test  results  with  model 
predictions.  The  test  data  are  limited  to  the  steering  performance  of 
a  selected  track -laying  vehicle  tested  at  three  different  soft-soil 
sites.  The  characteristics  of  the  soil  at  the  test  locations  are 
described  and  values  of  the  soil  model  material  constants  are  developed 
from  in  situ  direct  shear  measurements  in  Part  II.  The  test  procedure, 
prediction  methodology,  and  comparisons  of  model  predictions  with  test 
data  are  presented  in  Part  III.  Part  IV  contains  a  summary  and  recom¬ 
mendations  for  future  work.  In  Appendix  A  the  terrain-vehicle  model 
reported  by  Baladi  and  Rohani  (1979  and  1981)  is  extended  to  include 
the  treatment  of  sloping  terrain  under  nonuniform  (transient)  turning 
motion.  Description  and  use  of  a  direct  shear  device  for  measurement 
of  pertinent  soil  properties  are  documented  in  Appendix  B.  Soil 
classification  data  for  all  test  sites  are  given  in  Appendix  C. 


PART  II:  SITE  CHARACTERIZATION 


Background 

3.  One  of  the  most  important  engineering  properties  of  a  soil 
affecting  trafficability  is  its  in  situ  shear  strength.  The  shear 
strength  of  earth  materials  varies  greatly  for  different  types  of  soil 
and  is  dependent  on  the  confining  pressure  and  time  rate  of  loading 
(shearing).  This  dependence,  however,  varies  with  respect  to  the 
soil's  cohesive  and  frictional  properties.  It  has  been  found  experi¬ 
mentally  that  the  shear  strength  of  purely  cohesive  soils  (such  as  a 
saturated  plastic  clay)  is  relatively  independent  of  the  confining 
stress,  but  strongly  affected  by  the  time  rate  of  shearing.  On  the 
other  hand,  the  shear  strength  of  purely  frictional  soils  (such  as  a 
dry  clean  sand)  is  found  to  oe  relatively  independent  of  time  rate  of 
loading,  but  strongly  dependent  on  confining  pressure.  The  shearing 
resistance  of  most  soils,  however,  is  due  to  both  frictional  and 
cohesive  components.  An  appropriate  test  for  determining  shear  strength 
for  application  in  mobility  studies  is  a  direct  shear  test  conducted  in 
situ  on  the  soil  surface.  A  field  direct  shear  device  has  been  devel¬ 
oped  at  the  WES  for  such  applications  and  is  documented  in  Appendix  B. 
This  device  was  used  to  measure  the  in  situ  shear  strength  of  the  soil 
at  each  test  location. 


Test  Sites 


4.  The  mobility  tests  were  conducted  at  two  sites  (test  sites  7A 
and  7B)  on  a  floodplain  north  of  Redwood,  Mississippi,  and  on  a  hydrau¬ 
lically-filled  dredge  spoil  area  (test  site  8)  within  the  WES  reserva¬ 
tion  at  Vicksburg,  Mississippi  (Schreiner  and  Green,  1980).  The 
floodplain  soil  at  Redwood  is  a  soft,  plastic  clay  classified  as  CH 
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according  to  the  Unified  Soil  Classification  System*  (USCS) ;  however, 
test  site  7A  was  on  plowed  ground  while  test  site  7B  was  on  natural 
ground.  The  dredged  soil  at  test  site  8  is  a  lean,  brown  clay  classi¬ 
fied  as  CL  according  to  the  USCS.  Six  test  series  were  conducted  at 
these  sites  during  June  1979.  Table  1  lists  the  location  and  date  of 
each  test  series  and  also  includes  the  average  rating  cone  index  (RCI) , 
a  measure  of  the  soil's  remolded  resistance  to  penetration  by  the 
standard  WES  cone  penetrometer  (Smith,  1964),  for  each  test  series. 

Soil  classification  data  for  the  test  locations  are  given  in  Appendix  C. 
It  is  noted  from  Table  1  that  within  each  test  site  the  strength  of  the 
soil  varies  for  different  test  locations.  For  example,  test  series  107- 
111  and  139-144  were  both  conducted  at  test  site  7B,  but  the  strength  of 
the  soil  at  the  two  test  locations  is  different. 

Direct  Shear  Test  Results 


5.  Eighty-three  direct  shear  tests  were  conducted  at  the  test 
sites  using  the  procedure  outlined  in  Appendix  B.**  Both  slow  and  fast 
tests  were  conducted  to  ascertain  the  sensitivity  of  the  shear  strength 
of  the  material  to  the  rate  of  deformation.  In  the  case  of  fast  tests, 
the  soil  specimens  were  sheared  at  a  strain  rate  of  approximately  0.125 
to  0.25  per  sec.  This  range  of  strain  rates  is  compatible  with  the 
average  strain  rate  experienced  by  the  soil  particles  under  the  track  of 
the  test  vehicle  during  steering.  This  strain  rate  is  estimated  from 
slip  velocity  calculations  to  be  on  the  order  of  0.1  per  sec.  The  slip 
velocity  calculations  are  based  on  actual  measurements  of  track  velocity 
and  vehicle  speed  during  steering. 


*  The  Unified  Soil  Classification  System  is  described  in  Technical 
Memorandum  No.  3-357  by  the  U.  S.  Army  Engineer  Waterways  Experiment 
Station  (1953) . 

**  Raw  data  from  these  tests  are  available  at  the  WES.  Only  a  summary 
of  the  synthesized  data  necessary  for  model  predictions  is  presented 
in  this  report. 


6.  Direct  shear  tests  were  conducted  using  applied  normal  loads 
of  8.6,  36.6,  65.7,  and  122.7  lb*  corresponding,  respectively,  to 
normal  stresses  of  0.54,  2.29,  4.11,  and  7.67  psi.  The  ground  contact 
pressure  for  the  test  vehicle  is  5.71  psi.  Therefore,  the  range  of 
normal  stresses  used  in  the  test  program  is  applicable  to  the  test 
vehicle  of  interest.  Figure  1  shows  typical  load-deformation  data 
obtained  from  the  direct  shear  device.  As  indicated  in  Figure  1  in  the 
case  of  fast  tests,  the  soil  specimens  were  monotonically  sheared  to 
failure,  whereas  the  slow  tests  experienced  several  load-unload  cycles 
of  deformation.  The  load-deformation  data  were  used  to  construct 
representative  failure  curves  for  the  six  test  locations  indicated  in 
Table  1.  The  representative  failure  curves  are  portrayed  in  Figures 
2-7  as  plots  of  peak  shear  load  versus  normal  load  for  both  the  slow 
and  fast  tests.  For  load-deformation  curves  that  did  not  exhibit  a 
distinct  peak,  the  value  of  shear  load  at  15  percent  strain  (correspond¬ 
ing  to  0.6  in.  of  deformation)  was  selected  to  construct  the  failure 
curves.  The  data  indicate  that  the  increase  in  the  shear  load  beyond 
15  percent  strain  is  negligible.  The  average  values  of  wet  density  and 
water  content  given  for  each  of  the  test  locations  in  Figures  2-7 
correspond  to  the  average  of  the  wet  density  and  water  content  of  the 
soil  specimens  tested  at  each  location.  The  failure  curves  in  these 
figures  clearly  demonstrate  that  the  shear  strength  of  the  soil  for  all 
three  sites  is  sensitive  to  the  rate  of  deformation.  It  is  further 
observed  from  Figures  2-7  that  the  shear  strength  of  the  CL  soil  at 
site  8  is  considerably  higher  than  the  shear  strength  of  the  CH  soil  at 
sites  7A  and  7B.  The  same  relative  difference  in  shear  strength 
between  the  two  materials  is  also  reflected  in  the  RCI  readings  listed 
in  Table  1. 


Material  Constants  for  Soil  Model 


7.  As  pointed  out  previously,  the  purpose  of  the  direct  shear 


*  A  table  of  factors  for  converting  U.  S.  customary  units  of  measure¬ 
ments  to  metric  (SI)  units  is  presented  on  page  3. 
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tests  was  to  determine  the  in  situ  shear  strength  of  the  material  for 
site  characterization  within  the  framework  of  the  soil  model  described 
in  Appendix  A.  Such  a  characterization  is  necessary  in  order  to  use 
the  WES  terrain-vehicle  interaction  model  for  predicting  the  steering 
performance  of  the  vehicle  on  the  particular  terrain  of  interest.  The 
soil  model  contains  six  material  constants.  Three  of  the  material 
constants  (A,  M,  and  N)  describe  the  static  failure  envelope  of  the 
mater.' al  (see  Equation  A1  and  Figure  Al,  Appendix  A);  two  parameters  (A 
and  Cj)  define  the  contribution  to  cohesive  strength  of  the  material 
due  to  dynamic  loading  (see  Equation  A2  and  Figure  Al);  and  one  parame¬ 
ter  (G)  defines  the  initial  shear  stiffness  coefficient  of  the  soil 
(st--1  Equation  A3  and  Figure  A2) .  The  first  step  in  determining  the 
numerical  values  of  the  material  constants  A,  M,  and  N  is  to  convert 
the  slow  test  failure  :urves  in  Figures  2-7  to  shear  strength 
versus  normal  stress  a  relations  by  simply  dividing  the  shear  load 
and  norm;  j.oad  by  the  cross-sectional  area  of  the  soil  specimen  (4-  by 
4-in.  specimen).  The  second  step  involves  fitting  Equation  Al  to  the 
resulting  versus  a  curves.  Figures  8-13  portray  the  experimental 

versus  o  curves  and  the  corresponding  model  behavior  for  each  of 
the  test  locations.  It  is  observed  from  these  figures  that  the  agree¬ 
ment  between  field  measurements  and  model  behavior  is  excellent.  The 
parameter  C<j  ,  corresponding  to  the  increase  in  soil  cohesion  due  to 
dynamic  loading  (maximum  loading  rate  of  interest),  can  be  determined 
from  Figures  2-7.  Basically,  corresponds  to  the  difference  between 

the  fast  and  slow  failure  curves  at  zero  normal  load  divided  by  the 
crosr -sectional  area  of  the  specimen.  The  parameter  A  defines  the 
rate  of  increase  in  soil  cohesion  due  to  deformation  velocity  (Equation 
A2).  If  shear  test  data  were  available  for  several  rates  of  deforma¬ 
tion,  the  value  of  A  would  be  determined  by  fitting  Equation  A2  (for 
o  =  0)  to  a  plot  of  cohesion  versus  deformation  velocity,  as  indicated 
in  Figure  B5,  Appendix  B.  In  the  absence  of  such  information,  however, 

A  is  treated  as  a  fitting  parameter;  i.e.,  the  value  of  A  is  deter¬ 
mined  on  the  basis  of  fitting  Equation  A4  to  the  stress-deformation 
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curves  from  fast  tests.  Before  fitting  Equation  A4  to  the  stress- 
deformation  curves,  the  value  of  the  initial  shear  stiffness  coefficient 
G  must  be  known.  The  parameter  G  is,  in  general,  dependent  on  the 
normal  stress  and  the  deformation  velocity  (Figure  B5) .  Within  the 
framework  of  the  present  soil  model,  however,  G  is  assumed  to  be  a 
constant  and  independent  of  either  the  normal  stress  or  the  deformation 
velocity.  Parametric  studies  conducted  with  the  terrain-vehicle  model 
have  indicated  that  the  steering  performance  of  the  vehicle  is  only 
mildly  dependent  on  the  soil  parameter  G  .  Therefore,  it  is  not 
necessary  to  characterize  this  parameter  precisely.  The  procedure 
adopted  to  determine  G  is  first  to  plot  the  initial  slopes  of  the 
stress-deformation  curves  versus  normal  stress  from  both  the  slow  and 
fast  tests.  Then,  using  this  plot,  select  an  average  value  of  G  at  a 
normal  stress  level  corresponding  to  the  ground  contact  pressure  for 
the  vehicle  of  interest. 

8.  The  numerical  values  of  the  six  material  constants  were  deter¬ 
mined  for  all  the  test  locations  following  the  above  procedures  and  are 
summarized  in  Table  2,  To  demonstrate  the  validity  of  the  soil  model 
for  simulating  the  stress-deformation  response  of  the  material  from 
direct  shear  tests,  typical  experimental  data  from  fast  tests  are 
compared  in  Figures  14-19  with  the  corresponding  model  behavior.  The 
deformation  velocity  A  of  0.75  in. /sec  used  in  the  soil  model  calcu¬ 
lations  corresponds  to  an  average  velocity  for  the  fast  tests.  The 
comparisons  were  made  for  a  normal  stress  a  of  4.11  psi,  the  direct 
shear  test  value  that  was  closest  to  the  ground  contact  pressure  of  the 
test  vehicle  (5.71  psi).  Comparisons  of  the  field  data  with  model 
behavior  in  Figures  14-19  indicate  that  the  soil  model  is  capable  of 
simulating  the  pertinent  features  of  the  shear  stress-deformation 
response  of  the  material  very  accurately.  Some  of  the  field  measure¬ 
ments  (e.g.,  Figures  18  and  19)  exhibit  some  degree  of  strain-softening 
that  cannot,  be  simulated  with  the  current  version  of  the  soil  model. 
However,  lor  the  present  application  of  the  model,  such  strain-softening 
behavior  is  not  a  significant  phenomenon. 
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PART  III:  COMPARISON  OF  THEORETICAL  PREDICTIONS  WITH  TEST  RESULTS 


Test  Procedure 


9.  The  field  tests  consisted  or  a  number  of  circular-turn  tests; 
details  of  the  test  procedure  are  documented  in  a  study  by  Green  (in 
preparation).  The  principal  objectives  of  the  circular-turn  tests 
were:  (a)  determine  the  effects  of  turning  radius  on  turning  perform¬ 
ance  in  terms  of  vehicle  speed  and  power  requirements,  (b)  determine 
the  effects  of  soil  strength  on  turning  performance,  and  (c)  develop  a 
data  base  to  check  the  accuracy  and  range  of  application  of  the  WES 
terrain-vehicle  interaction  model.  Basically,  each  test  involved 
running  the  tracked  vehicle  in  a  circular  path  by  first  accelerating 
the  vehicle  to  a  maximum  speed  (controlled  by  either  the  available 
power  or  the  actual  physical  stability  of  the  vehicle)  and  then  con¬ 
tinuously  turning  it  in  a  near  steady-state  condition.  The  tracked 
vehicle  used  for  the  field  tests  is  an  armored  personnel  vehicle  with 
the  characteristics  listed  in  Table  3.  The  actual  data  collected 
during  each  test  consisted  of  time  histories  of  (a)  sprocket  rpm  and 
torque,  (b)  turning  radius,  and  (c)  velocity  of  the  vehicle.  Then  the 
track  velocities,  lateral  acceleration,  and  power  consumption  were 
calculated  for  each  test  using  these  measurements  and  appropriate 
equations  (Figure  20).  In  addition  to  the  above  data,  stopwatch  times 
for  each  revolution  and  posttest  measurements  from  the  center  of  the 
circle  to  the  inner  and  outer  track  ruts  were  obtained  to  calculate  an 
average  effective  vehicle  speed  and  turning  radius  for  each  circular 
path.  This  information  provided  a  check  on  the  turning  instrumentation. 
Also,  in  order  to  determine  the  coefficient  of  rolling  resistance  (( 

(see  Equation  A46) ,  acceleration/deceleration  (AC/DC)  tests  were  con¬ 
ducted  at  each  test  location  (Green  (in  preparation)).  Table  4  lists 
the  resulting  values  of  for  the  six  test  locations. 
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Prediction  Methodology 

10.  The  prediction  methodology  associated  with  the  WES  terrain- 
vehicle  interaction  model  is  illustrated  in  Figure  21.  The  model  input 
consists  of  three  separate  sets  of  data.  The  first  set  of  data 
describes  the  mechanical  properties  of  the  terrain  within  the  framework 
of  the  soil  model  described  in  Appendix  A.  Table  2  presents  the  data 
necessary  for  this  set  of  input  parameters  for  all  test  locations.  The 
second  set  of  data  describes  the  characteristics  of  the  vehicle  and  the 
coefficient  of  rolling  •  esistance  for  each  test  location.  Tables  3  and 
b  outline  such  information  for  the  vehicle  of  interest  for  all  the  test 
locations.  The  third  set  of  input  data  describes  the  conditions  by 
which  one  would  drive  the  model.  There  are  several  combinations  of 
driving  conditions  that  can  be  used  depending  on  the  nature  of  the 
particular  problem  at  hand  (Baladi  and  Rohani,  1979).  For  the  particu¬ 
lar  application  of  interest  in  this  report,  the  driving  conditions 
consist  of  the  time  histories  of  the  inner  and  outer  track  velocities, 
which  are  obtained  from  actual  field  measurements.  For  these  specified 
driving  conditions,  the  model  outputs  the  time  histories  of  vehicle 
velocity,  slip  velocities,  lateral  acceleration,  power  at  sprockets, 
turning  radius,  and  offset.  These  quantities  can  then  be  compared  with 
corresponding  field  data  in  order  to  determine  the  accuracy  of  the 
model.  As  indicated  in  Figure  20,  the  field  data  for  the  circular-turn 
tests  consist  of  vehicle  velocity,  turning  radius,  lateral  acceleration, 
and  power  at  the  sprockets.  The  average  slip  velocity  can  also  be 
calculated  from  the  direct  field  measurements  if  desired. 

11.  In  using  the  time  histories  of  track  velocities  to  drive  the 
terrain-vehicle  interaction  model,  the  field  measurements  were  filtered 
and  fitted  by  analytical  expressions  because  the  field  data  included 
artificial  high-frequency  oscillations  that  are  believed  to  be  instru¬ 
mentation  "noise."  These  oscillations  lead  to  frequent  crossings  of 
the  two  track  velocity-time  histories,  which  are  unrealistic  for  the 
circular-turn  tests  of  interest.  To  filter  the  data,  each  track 
velocity-time  history  was  first  numerically  integrated.  The  integrated 
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results  were  plotted  versus  time  and  were  represented  by  a  Fourier 
series  containing  21  terms.  The  Fourier  series  was  then  differentiated 
to  yield  an  expression  for  the  track  velocity.  The  differentiated 
analytical  expression  for  the  track  velocity  has  the  following  form 


=  Ai  + 


Bl  + 


21 
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i=2 


27r(i-l) 
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27T(j-l)t 


(1) 


where 

v^  =  track  velocity 

T  =  total  time  indicating  the  duration  of  the  event 
t  =  time 

A  ,  (i  =  1,...,21)  =  constants 

The  numerical  values  of  the  coefficients  A  and  are  given  in 

Tables  5-13  for  use  of  Equation  1  by  interested  individuals.  Compari¬ 
sons  of  the  filtered  data  and  field  measurements  for  all  the  tests  are 
presented  in  the  next  section. 

Theoretical  Predictions 


12.  The  experimental  program  consisted  of  35  circular  tests 
divided  into  six  test  series  (Table  1).  Unfortunately,  due  to  instru¬ 
mentation  problems,  complete  data  (as  described  in  Figure  20)  were 
collected  for  only  16  tests.  These  tests  are  referred  to  as  "good 
tests."  The  data  for  other  tests  were  either  partially  complete  or  of 
poor  quality.  The  theoretical  predictions  reported  in  this  section 
include  the  16  good  tests  (all  from  sites  7A  and  7B)  and  one  test  with 
poor  quality  data  (from  site  8).  The  16  tests  consist  of  tests  99,  100, 
101,  102,  107,  108,  109,  110,  111,  117,  118,  133,  134,  139,  140,  and 
141.  Test  157  (from  test  series  157-159)  was  selected  as  an  example  of 
a  test  with  poor  quality  data.  Two  types  of  predictions  are  made  for 
the  tests  above  and  are  discussed  in  the  following  paragraphs. 
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Complete  transient  simulation 

13.  To  demonstrate  the  ability  of  the  terrain-vehicle  model  for 
predicting  transient  motion,  the  response  of  the  vehicle  was  simulated 
for  the  entire  test  event.  The  results  of  these  simulations  are 
portrayed  in  five  figures  for  each  test.  The  first  figure  consists  of 
the  time  histories  of  the  outer  and  inner  track  velocities  that  con¬ 
stitute  the  driving  conditions  for  the  model.  This  figure  includes 
both  the  field  measurements  and  the  filtered  data  (see  paragraph  11), 
which  are  used  as  input  to  the  model.  The  next  three  figures  show  the 
time  histories  of  the  vehicle  speed,  lateral  acceleration,  and  power 
requirement,  and  include  both  the  field  measurements  and  the  correspond¬ 
ing  model  predictions.  The  fifth  figure  shows  the  trajectory  of  the 
center  of  gravity  of  the  vehicle  and  includes  both  field  measurements 
and  predictions.  Figures  22-106  show  the  simulation  results  for  the  17 
tests  indicated  in  paragraph  12.  The  figure  numbers  corresponding  to 
each  test  are  listed  in  the  following  tabulation: 


Figure  No, 


Test 

Track 

Vehicle 

Lateral 

Power 

No. 

Velocity 

Velocity 

Acceleration 

Requirement 

Trajectory 

99 

22 

23 

24 

25 

26 

100 

27 

28 

29 

30 

31 

101 

32 

33 

34 

35 

36 

102 

37 

38 

39 

40 

41 

107 

42 

43 

44 

45 

46 

108 

47 

48 

49 

50 

51 

109 

52 

53 

54 

55 

56 

110 

57 

58 

59 

60 

61 

111 

62 

63 

64 

65 

66 

117 

67 

68 

69 

70 

71 

3.18 

72 

73 

74 

75 

76 

133 

77 

78 

79 

80 

81 

134 

82 

83 

84 

85 

86 

139 

87 

88 

89 

90 

91 

140 

92 

93 

94 

95 

96 

141 

97 

98 

99 

100 

101 

157* 

102 

103 

104 

105 

106 

*  Representative  test  with  poor  quality  data. 


14 


The  measured  time  histories  basically  manifest  two  types  of  oscilla¬ 
tions — low  frequency  and  high  frequency.  The  high-frequency  oscilla¬ 
tions,  as  pointed  out  earlier,  are  primarily  due  to  instrumentation, 
and  no  physical  interpretation  should  be  attached  to  them.  The  low- 
frequency  oscillations,  however,  are  real  and  are  partly  due  to  driver 
response  (i.e.,  sudden  acceleration  and  deceleration  of  the  vehicle  due 
to  steering  correction  to  maintain  a  prescribed  circular  path)  and 
partly  due  to  surface  roughness,  inhoinogeneity  of  surface  materials, 
and  the  fact  that  the  ground  is  not  an  ideally  flat,  level  surface. 
Because  of  these  factors,  it  is  not  possible  to  maintain  the  vehicle  in 
a  perfect  steady-state  mode  of  motion.  Therefore,  in  comparing  the 
field  data  with  model  predictions,  it  is  the  overall  response  that  must 
be  considered,  not  the  peaks  and  valleys  of  the  oscillatory  records. 
With  this  in  mind,  the  comparisons  between  model  predictions  and  field 
measurements  are  very  reasonable  for  the  "good  tests"  (Figures  22-101) . 
For  several  of  these  tests,  however,  the  predicted  power  requirement  is 
lower  than  the  field  data  (e.g.,  Figures  30  and  35).  Examination  of 
the  field  notes  indicated  that  there  was  excessive  mud  buildup  inside 
the  tracks  of  the  vehicle  for  these  tests.  The  mud  buildup  (which  is 
not  simulated  by  the  model)  would  require  additional  power  to  steer  the 
vehicle.  The  results  of  test  157  (poor  quality  data)  are  presented  in 
Figures  102-106.  It  is  obvious  from  these  figures  that  except  for 
turning  radius  the  field  measurements  are  of  poor  quality. 

Steady-state  simulation 

1A.  As  pointed  out  in  the  previous  paragraph,  it  was  not  possible 
to  maintain  the  vehicle  in  a  perfect  steady-state  mode  of  motion  during 
the  entire  test  event.  However,  it  is  possible  to  select  a  small  time 
window  for  each  test  where  the  motion  of  the  vehicle  can  be  reasonably 
approximated  as  steady  state.  The  steady-state  version  of  the  terrain- 
vehicle  interaction  model  (Appendix  A)  can  then  be  used  to  simulate 
such  motions.  Such  simulations  were  documented  for  tests  107-111 
(Baladi  and  Rohani,  1981);  for  ease  of  reference,  the  results  are  again 
presented  in  this  report.  The  time  windows  and  the  corresponding  test 
data  consisting  of  the  steering  ratio  e  ,  vehicle  velocity  v  , 
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2 

turning  radius  Rq  ,  laterai  acceleration  v  /RQg  ,  and  total  power  PT 
for  conducting  the  steady-state  simulations  are  given  in  the  following 
tabulation  for  all  of  these  tests: 


Test 

No. 

Time 

sec 

Steering 

Ratio 

e 

Vehicle 
Speed 
v,  mph 

Turning 

Radius 

R  ,  ft 

0 

Lateral 

Acceleration 

v/,/R0g 

Power 
PT,  hr 

107 

60 

66 

1.08 

24.88 

155.77 

0.265 

210 

108 

81 

84 

1.29 

14.93 

63.69 

0.234 

199 

1  AQ 

37 

42 

1.12 

15.73 

103.86 

0.16 

177 

iuy 

68 

71 

1.10 

18.22 

121.26 

0.18 

184 

110 

55 

59 

1.26 

16.08 

83.40 

0.21 

200 

111 

48 

51 

1.57 

11.79 

33.74 

0.275 

211 

Comparisons  of  the  model  predictions  with  experimental  data  are 
presented  in  Figures  107-112  for  turning  radius  versus  steering  ratio, 
inner  track  velocity  versus  turning  radius,  outer  track  velocity  versus 
turning  radius,  vehicle  speed  versus  turning  radius,  power  requirement 
versus  turning  radius,  and  lateral  acceleration  versus  turning  radius, 
respectively.  The  model  predictions  in  Figures  108-112  are  based  on 
both  power  cutoff  and  preliminary  stability  criteria  (see  Appendix  A) . 
The  turning  radius-steering  ratio  relation  shown  in  Figure  107,  however, 
is  unique  for  a  given  vehicle  and  soil  condition.  The  power  cutoff,  as 
indicated,  is  controlled  by  the  available  power.  The  preliminary 
stability  criteria  are  based  on  (Baladi  and  Rohani,  1979): 

a.  Rapid  change  in  the  slip  velocity  of  the  inner  or  the  outer 
track. 

_b.  The  pivot  point  falling  outside  the  front  edge  of  the 

track-ground  contact  area  (i.e.,  the  offset  equals  0.5  L 
when  the  center  of  gravity  and  center  of  geometry  of  the 
vehicle  coincide) . 

£.  Rapid  decrease  or  increase  in  the  turning  radius. 

These  stability  conditions  usually  take  place  at  different  vehicle 
velocities.  The  unstable  vehicle  velocity  is  chosen  as  the  minimum 
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velocity  at  which  any  of  these  conditions  occur.  For  comparison  with 
~he  experimental  data,  the  lower  vehicle  velocity  corresponding  to 
either  the  stability  criteria  or  the  power  cutoff  condition  must  be 
selected.  As  indicated  in  Figure  110,  for  the  turning  radii  of  34  ft 
(test  111)  and  64  ft  (test  108),  stability  criteria  control  the  velocity 
of  the  vehicle.  For  the  turning  radii  of  83  ft  (test  110),  104  and  121 
ft  (test  109),  and  156  ft  (test  107),  the  velocity  of  the  vehicle  is 
controlled  by  the  available  power.  With  this  in  mind,  the  experimental 
data  in  Figures  107-112  compare  very  favorably  with  the  corresponding 
model  predictions.  This  fact  is  particularly  true  in  the  case  of  track 
velocities  and  vehicle  speed  (Figures  108-110).  Slight  observable 
differences  between  the  data  and  model  predictions  in  Figures  107-112 
should  be  expected  because  of  the  small  deviations  in  the  test  condi¬ 
tions  from  the  steady-state  mode  of  motion. 


\ 
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PART  IV:  SUMMARY  AND  RECOMMENDATIONS 

15.  In  1978,  the  WES  developed  a  mathematical  model  to  predict  the 
steering  performance  of  tracked  vehicles.  In  1979,  field  tests  of  a 
selected  track-laying  vehicle  were  conducted  at  plowed  and  unplowed 
sites  on  a  floodplain  north  of  Redwood,  Mississippi,  and  on  a  dredge 
spoil  area  within  the  WES  reservation  at  Vicksburg,  Mississippi.  In 
this  report,  the  test  measurements  are  compared  with  calculated  results 
in  order  to  validate  the  WES  terrain-vehicle  interaction  model  for 
predicting  tracked-vehicle  performance  on  soft  soils. 

16.  The  floodplain  soil  at  Redwood  is  a  soft,  plastic  clay  (CH) ; 
the  dredged  soil  at  the  WES  is  a  lean,  brown  clay  (CL).  Eighty-three 
direct  shear  tests  were  conducted  at  the  test  sites.  The  shear  strength 
of  the  soils  at  all  three  sites  was  sensitive  to  the  rate  of  deforma¬ 
tion;  therefore  a  soil  model  with  a  rate-dependent  nonlinear  failure 
envelope  was  developed.  The  soil  model  contains  six  material  constants; 
values  for  these  constants  were  defined  for  each  site  by  fitting  the 
model  to  the  in  situ  test  results.  The  agreement  between  the  soil  model 
fits  and  the  field  data  is  excellent. 

17.  A  total  of  35  circular-turn  tests  were  conducted  at  the  three 
sites.  Due  to  instrumentation  problems,  however,  only  16  were  con¬ 
sidered  "good  tests"  for  validating  the  tracked-vehicle  model.  Theoret¬ 
ical  predictions  were  made  for  these  16  "good  tests,"  all  of  which  were 
conducted  at  the  two  sites  on  the  floodplain  north  of  Redwood.  Theoret¬ 
ical  predictions  were  also  made  for  one  test  with  poor  quality  data 
from  the  dredged  soil  site.  Two  types  of  predictions  were  made.  The 
first  type  was  made  using  the  transient  version  of  ".he  model;  the 
second  was  made  using  the  steady-state  version  of  the  model. 

18.  The  input  driving  conditions  for  the  16  transient  predictions 
consisted  of  time  histories  of  inner  and  outer  track  velocities 
obtained  from  actual  field  measurements.  To  eliminate  artificial  high- 
frequency  oscillations  and  frequent  unrealistic  crossings  of  the  two 
track  velocity-time  histories,  the  measured  track  velocities  were 
filtered  using  Fourier  analysis.  The  model  output  consisted  of  time 
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histories  of  vehicle  speed,  lateral  acceleration,  power  required  at  the 
sprockets,  and  the  trajectory  of  the  center  of  gravity  of  the  vehicle. 
Comparisons  with  the  corresponding  experimental  results  are  very 
favoraale,  indicating  that  the  model  is  capable  of  predicting  the 
steering  performance  of  tracked  vehicles  on  soft  soil. 

19.  Although  the  vehicle  tests  were  intended  to  be  constant 
velocity  or  steady-state  turn  tests,  the  measured  track  velocities  were 
not  constant  partly  because  of  (a)  driver  response  (i.e.,  sudden  ac¬ 
celeration  or  deceleration  of  the  vehicle  due  to  steering  correction  to 
maintain  a  prescribed  circular  path),  and  (b)  terrain  roughness,  in¬ 
homogeneity  of  surface  materials,  pj.us  the  fact  that  the  terrain  is 

not  an  ideally  flat,  level  surface.  Therefore,  even  "steady-state" 
turn  tests  should  be  simulated  using  the  transient  version  of  the 
model. 

20.  While  it  was  not  possible  to  maintain  the  vehicle  in  a  perfect 
steady-state  mode  of  motion  (i.e.,  constant  track  velocities)  during  the 
entire  test  event,  steady-state  conditions  did  exist  for  short  periods 
of  time  during  each  test.  Steady-state  calculations  were  performed  for 
selected  time  windows  in  five  of  the  "good  tests."  The  input  consisted 
of  a  mean  radius  for  each  test;  the  output  consisted  of  a  maximum 
vehicle  velocity  based  on  either  specified  stability  criteria  or  power 
available  at  sprockets.  These  predictions  also  correlated  very  favorably 
with  the  corresponding  test  data. 

21.  To  further  validate  the  accuracy  and  range  of  application  of 
the  model,  data  are  needed  from  (a)  maneuvering  tests  using  a  number  of 
tracked  vehicles  in  which  the  vehicle  cnaracteristics  (such  as  weight, 
track  length,  and  tread)  are  varied,  (b)  tests  conducted  on  sloping 
terrains,  and  (c)  tests  conducted  on  soils  other  than  a  soft,  plastic 
clay  (e.g.,  purely  cohesionless  soils  such  as  dry  sand  and  very  hard 
surfaces  such  as  pavements) . 
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Table  1 

Test  Site,  Soil  Type,  and 
Soil  Strength  for  Each  Test  Series 


Test 

Test* 

Soil 

Average 

Soil  Strength 

RCI 

Series 

Site 

Type 

0-6  in. 

6-12  in. 

Date 

98-102 

7A 

CH 

40 

69 

1 

June 

1979 

107-111 

7B 

CH 

24 

51 

5 

June 

1979 

117-122 

7A 

CH 

57 

93 

7 

June 

1979 

131-134 

7A 

CH 

61 

96 

8 

June 

1979 

139-144 

7B 

CH 

54 

85 

14 

June 

1979 

157-159 

8 

CL 

210 

198 

28 

June 

1979 

*  Site  7  A  is  a  plowed  ground;  site  7B  is  natural  ground; 
and  site  8  is  a  hydraulic  fill. 


Table  2 

Values  of  Soil  Model  Material  Constants  for 
Each  Test  Series 


Test 

Series 

•V 

psi 

M 

psi 

N 

1/psi 

Cd 

psi 

/; 

sec/in. 

G 

psi/in. 

98-102 

3.5 

1.94 

0.3 

1.25 

10.0 

150.0 

107-111 

5.0 

3.75 

0.22 

1.10 

10.0 

200.0 

117-122 

7.8 

6.36 

0.1 

1.75 

10.0 

200.0 

131-134 

5.4 

3.65 

0.23 

2.14 

10.0 

125.0 

139-144 

6.5 

4.94 

0.09 

1.75 

10.0 

270.0 

157-159 

24.0 

19.94 

0.02 

2.64 

10.0 

250.0 
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Table  3 

Characteristics  of  Vehicle  Used  for  Turn  Tests 


Weight  (W) 

= 

18,000 

Ibm 

Track  Length  (L) 

= 

105 

in. 

Track  Width  (D) 

= 

15 

in. 

Tread  (B) 

= 

90 

in. 

Height  of  the  center  of  gravity  (H) 

= 

35.7 

in. 

Location  of  the  center  of  gravity  measured  from 
the  geometrical  center  of  the  vehicle  (C^) 

= 

•  0 

in. 

Distance  between  two  adjacent  wheels  (l) 

= 

26.25 

in. 

Mass  moment  of  inertia  (Iz) 

= 

92,000 

...  2 
lbf-in.-sec 

Approach  angle  (8a) 

= 

30 

deg 

Departure  angle  (0^) 

= 

30 

deg 

\ 


Table  4 

Rolling  Resistance  Measured  from  AC/DC  Tests 


Coefficient  of  Rolling  Resistance 


Test  Series 

98-102 

0.19 

107-111 

0.20 

117-122 

0.16 

131-134 

0.14 

.139-144 

0.16 

157-159 

0.125 
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Table  5 

Numerical  Values  of  the  Coefficients  and  B* 


Equation  1)  for  Tests  99  and  100 


TEST  NUMBER  99 


TEST  NUMBER  100 


OUTER  TRACK 


INNER  TRACK 


OUTER  TRACK 


INNER  TRACK 


13.578 

69.935 

18.701 

4.137 

6.782 


7  1.432 

8  1.704 
0.807 

10  0.735 

11  0.980 

12  0.889 

13  0.675 

14  0.251 

15  0.377 

16  0.274 

17  0.408 

18  0.245 

19  0.563 

20  0.407 

21  0.248 


9.2699 
46.269 
14.431 
2.738 
5.189 
1.300 
1.120 
1.347 
0.626 
0.416 
0.823 
0.491 
0.131 
0.649 
.155  0.492 

.181  0.284 

.067  0.378 

.050  0.231 

.050  0.372 

.311  0.013 

0.149 


10.543  0 

59.974  -37.856 
6.788  -15.614 
5.167  -6.928 

1.060  -6.089 

2.457  -2.182 

-0.073  -3.239 

0.270  -1.512 

-0.847  -0.421 

0.151  -1.037 

< 

-0.544  -0.120 


-0.330 


0.078 


-0.354 


-0.468  0.010 

-0.232  -0.142 

0.246  0.018 

-0.119  0.073 

-0.106  0.144 


0.103 


-0.060  -0.080 
0.154  -0.077 

0.081  -0.123 


-25.739 

-  8.648 

-  3.914 

-  3.669 

-  1.522 

-  1.725 

-  1.013 

-  0.492 

-  0.783 

-  0.014 

-  0.427 

-  0.189 

-  0.148 

-  0.045 
0.181 

-  0.010 

-  0.145 
0.103 

-  0.023 

-  0.032 
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Table  6 

Numerical  Values  of  the  Coefficients  and 
(Equation  1)  for  Tests  101  and  102 


TEST  NUMBER  101  1 

l 

t 

s 

|  TEST  NUMBER  102 

OUTER  TRACK 

INNER  TRACK  | 

|  OUTER  TRACK 

INNER  TRACK 

mm 

n 

mm 

13.619 

0 

61.520 

-62.582 

50.982 

17.260 

4.926 

15.341 

7.404 

-  6.457 

6.167 

5.600 

-  6.986 

4.318 

1.030 

-  0.650 

1.255 

1.387 

-  2.295 

0.994 

1.717 

0.200 

1.478 

1.354 

0.086 

1.201 

1.516 

-  1.322 

0.920 

0.853 

-  1.094 

0.832 

0.655 

-  0.267 

0.385 

0.324 

-  0.339 

0.329 

0.579 

-  0.043 

0.567 

0.702 

-  0.255 

0.361 

0.220 

-  0.509 

0.063 

0.151 

-  0.405 

0.115 

0.115 

-  0.179 

0.074 

0.045 

0.174 

0.063 

-  0.080 

0.347 

0.207 

0.008 

0.222 

-51.107 

4.033  I  41.238 
-  5.522 


-0.270 

-1.951 

0.160 

-0.040 

-1.055 

-0.685 

-0.237 


-0.066 

-0.413 


-0.385 

0.069 


-0.020 


-0.130 


Bi 

HI 

HI 

0 

11.961 

0 

-84.309 

48.839 

-76.819 

-  7.026 

38.371 

-  6.153 

-  3.247 

0.497 

-  3.042 

0.988 

2.938 

0.989 

4.736 

-  2.617 

-  0.446 

4.196 

-  0.230 

-  1.386 

-0.328 

-  1.412 

mm 

1.571 

-  0.097 

-  0.573 

1.936 

-  0.479 

-  0.365 

0.693 

-  0.297 

-  0.410 

0. 59i,'  1 

1  -  0.172 

-  0.561 

0.57] 

-  0.490 

0.028 

0.852 

0.002 

-  0.169 

0.569 

-  0.141 

-  0.254 

0.431 

-  0.228 

-  0.326 

0.368 

-  0.340 

-  0.393 

0.425 

-  0.399 

-  0.130 

0.375 

-  0.053 

-  0.083 

0.208 

-  0.149 

-  0.264 

0.162 

-  0.160 

Table  7 

Numerical  Values  of  the  Coeff  icients_Ajand_Bi 
- -^nation  1  )  f*r  Tests  107  and  108 


Table  8 


Numerical  Values  of  the  Coefficients  Aj  and 
(Equation  1  )  for  Tests  109  and  110 


I 

N 

TEST  NUMBER  109 

|  TEST  NUMBER  110 

H 

OUTER  TRACK 

INNER 

TRACK 

1 

OUTER  TRACK 

INNER  TRACK 

H 

D 

H 

■mm 

Bi 

BI 

B 

i 

14.031 

0 

12.784 

14.797 

0 

12.690 

0 

2 

52.868 

-70.877 

47.416 

56.488 

-56.478 

43.883 

-45.917 

3 

20.035 

-12.873 

17.624 

j  21.489 

-  4.846 

18.416 

-  5.903 

4 

7.470 

-  3.498 

7.394 

-  3.761 

6.959 

-  2.340 

5.494 

-  3.293 

5 

3.473 

-  6.441 

2.413 

-  5.862 

4.732 

-  1.224 

4.645 

-  1.795 

6 

1.486 

-  3.986 

1.340 

-  3.896 

j  4.012 

-  3.666 

1.470 

-  3.381 

7 

0.388 

-  2.978 

0.254 

-  2.572 

0.734 

-  1.285 

0.864 

-  0.676 

8 

0.277 

-  2.222 

0.061 

-  1.891 

!  1.061 

-  1.524 

0.710 

-  1.479 

9 

-0.210 

-  1.441 

-0.114 

-  1.113 

1.140 

-  1.151 

0.464 

-  1.028 

10 

-0.352 

-  1.389 

-0.617 

BBS 

-  1.469 

-0.226 

-  0.775 

11 

-0.263 

-  1.255 

-0.476 

-  1.186 

0.296 

0.046 

0.561 

-  0.000 

12 

-0.411 

-  0.543 

-0.492 

-  0.451 

0.493 

-  0.345 

0.409 

-  0.458 

13 

-0.197 

-  0.298 

-0.034 

-  0.378 

0.358 

-  0.200 

C.  244 

-  0.297 

14 

-0.333 

-  0.603 

-0.326 

-  0.448 

-  0.314 

0.083 

-  0.113 

15 

0.026 

-  0.425 

-0.139 

-  0.436 

-  0.069 

0.305 

-  0.281 

16 

-0.391 

-  0.496 

-0.373 

-  0.182 

0.243 

-  0.233 

0.073 

-  0.341 

17 

-0.050 

-  0.048 

-0.028 

-  0.272 

0.153 

-  0.264 

0.062 

-  0.235 

18 

-0.236 

-  0.118 

-0.089 

-  0.079 

0.185 

-  0.300 

0.079 

-  0.216 

19 

0.121 

-  0.234 

-0.135 

-  0.278 

0.141 

0.017 

0.154 

-  0.144 

20 

-0.152 

-  0.285 

-0.099 

-  0.027 

0.233 

-  0.397 

-0.090 

-  0.258 

21 

-0.076 

-  0.092 

-0.003 

-  0.069 

0.006 

-  0.227 

-0.012 

-  0.087 

27 


Table  9 


Numerical  Values  of  the  Coefficients  Aj  and  Bj 
(Equation  1  )  for  Tests  111  and  117 


TEST  NUMBER  111 


TEST  NUMBER  117 


OUTER 

TRACK 

INNER  TRACK 

|  OUTER 

TRACK 

INNER 

TRACK 

mm 

Bi 

mm 

Bi 

n 

n 

B 

B 

12.484 

0 

8.201 

0 

18.062 

0 

16.643 

0 

27.350 

-33.928 

13.408 

-21.948 

124.848 

-134.382 

111.637 

-119.956. 

8.983 

-  3.087 

5.617 

-  2.509 

38.878 

2.946 

35.272 

2.856 

2.492 

-  3.579 

1.134 

-  2.442 

4.516 

1.541 

5.043 

1.924 

2.235 

-  3.012 

1.840 

-  2.127 

1.554 

-4.525 

2.061 

-2.873 

1.641 

-  1.656 

1.420 

-  0.836 

11.284 

10.782 

10.923 

8.031 

0.835 

-  1.092 

0.316 

-  1.023 

3.952 

-5.554 

2.834 

-4.562 

0.727 

-  0.535 

0.604 

-  0.374 

4.976 

-1.941 

4.349 

-1.730 

0.365 

-  0.568 

0.319 

-  0.277 

1.823 

0.120 

1.789 

0.213 

0.231 

-  0.619 

0.093 

-  0.341 

0.452 

-0.366 

0.608 

0.122 

0.306 

-  0.321 

0.134 

-  0.010 

3.198 

0.596 

3.228 

-0.225 

0.120 

0.110 

0.315 

0.063 

1.50S 

-1.189 

1.027 

-0.971 

0.577 

-  0.147 

0.223 

-  0.220 

1.353 

-0.482 

1.170 

-0.464 

0.386 

0.001 

0.266 

-  0.115 

0.888 

-0.279 

0.874 

-0.187 

0.250 

-  Q.118 

0.074 

-  0.165 

0.882 

-0.340 

0.872 

-0.140 

0.323 

-  0.022 

0.184 

-  0.149 

0.847 

-0.157 

0.970 

-0.443 

0.443 

-  0.173 

0.155 

-  0.261 

0.863 

-0.404 

0.526 

-0.321 

0.032 

-  0.364 

-0.098 

-  0.068 

0.549 

-0.340 

0.421 

-0.306 

0.070 

-  0.037 

0.042 

-  0.055 

0.369 

-0.238 

0.323 

-0.173 

0.033 

-  Q. 134 

0.026 

-  0.108 

0.412 

-0.083 

0.377 

-0.046 

0.013 

-  0.104 

0.007 

-  0.004 

0.345 

-0.098 

0.488 

-0.221 

Table  10 

Numerical  Values  of  the  Coefficients  Ai  and  Bi 


(Equation  1)  for  Tests  118  and  133) 


Table  11 


Numerical  Values  of  the  Coefficients  Ai  and  Bi 
(Equation  1)  for  Tests  134  and  139 


I 

N 

D 

E 

Y 

TEST  NUMBER  134 

TEST  NUMBER  139 

OUTER  TRACK 

INNER  TRACK 

OUTER 

TRACK 

INNER 

TRACK 

A. 

B, 

A. 

B. 

A. 

B , 

A. 

B. 

1 

i 

1 

l 

l 

i 

X 

1 

1 

11.273 

0 

8.009 

0 

21.834 

0 

20.196 

0 

2 

19.954 

-28.321 

9.413 

-17.684 

122.426 

-95.596 

111.791 

89.443 

3 

8.482 

-  2.735 

4.809 

-  2.295 

28.839 

-  2.653 

26.801 

-2.998 

4 

2.604 

-  2.526 

1.361 

-  1.539 

13.627 

-  8.559 

12.271 

-7.919 

5 

2.435 

-  1.710 

1.246 

-  0.831 

4.825 

-  4.434 

4.396 

-3.921 

6 

1.311 

-  0.589 

1.097 

-  0.203 

3.692 

-  1.132 

3.395 

-1.094 

7 

1.060 

-  0.709 

0.677 

-  0.494 

2.761 

-  2.631 

2.484 

-2.326 

8 

0.854 

-  0.600 

0.552 

-  0.247 

1.110 

-  1.663 

1.067 

-1.418 

9 

0.700 

-  0.341 

0.410 

-  0.166 

1.111 

-  0.398 

1.138 

-0.391 

10 

0.284 

-  0.379 

0.371 

-  0.163 

1.470 

-  0.906 

1.344 

-0.904 

11 

0.280 

-  0.166 

0.163 

-  0.220 

0.676 

-  0.773 

0.659 

-0.816 

12 

0.167 

-  0.281 

0.068 

-  0.129 

0.795 

-  0.591 

0.687 

-0.647 

13 

0.267 

-  0.114 

0.175 

-  0.209 

0.411 

-  0.422 

0.403 

-0.347 

14 

0.066 

-  0.221 

0.001 

0.025 

0.393 

-  0.286 

0.414 

-0.345 

15 

0.050 

-  0.041 

0.169 

-  0.014 

0.411 

-  0.264 

0.347 

-0.264 

16 

0.157 

-  0.134 

0.006 

-  0.038 

0.471 

-  0.380 

0.327 

-0.409 

17 

0.084 

-  0.100 

0.081 

-  0.042 

0.176 

-  0.326 

0.099 

-0.194 

18 

0.009 

-  0.156 

0.086 

0.028 

0.241 

-  0.224 

0.171 

-0.149 

19 

0.053 

0.002 

0.054 

-  0.107 

0.259 

-  0.274 

0.211 

-0.236 

20 

0.092 

-  0.123 

0.008 

0.003 

0.234 

-  0.226 

0.153 

-0.215 

21 

0.004 

-  0.092 

0.024 

-  0.025 

0.091 

-  0.228 

0.080 

-0.114 
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Table  12 

Numerical  Values  of  the  Coefficients  and 
(Equation  1)  for  Tests  140  and  141 


TEST  NUMBER  140 


TEST  NUMBER  141 


OUTER  TRACK 


INNER  TRACK 


OUTER  TRACK 


INNER  TRACK 


Ai 

B. 

l 

A. 

l 

B. 

l 

Ai 

Bi 

A. 

l 

19.704 

0 

17.737 

0 

13.798 

o  ! 

11.291 

74.107 

-54.121 

64.482 

-48.583 

48.173 

-90.447 

36.294 

17.046 

-  0.236 

15.305 

-  0.309 

9.322 

-12.514 

7.749 

7.022 

-  3.780 

6.295 

-  3.159 

-  4.608 

-  8.309 

-3.391 

4.384 

-  2.138 

4.000 

-  1.987 

-  4.030 

-  4.228 

-3.187 

1.940 

-  1.793 

1.670 

-  1.649 

-  2.995 

-  3.521 

-2.532 

2.320 

-  1.831 

1.893 

-  1.758 

! 

-  1.108 

-  1.814 

-0.990 

0.626 

-  1.510 

0.439 

-  1.028 

0.117 

0.030 

0.174 

0.880 

-  0.209 

1.073 

-  0.230 

0.135 

0.417 

0.288 

0.866 

-  0.739 

0.701 

-  0.780 

0.231 

0.272 

0.196 

0.480 

-  0.637 

0.349 

-  0.584 

0.395 

-  0.323 

0.106 

0.409 

-  0.527 

0.298 

-  0.528 

-  0.170 

-  0.707 

-0,339 

0.352 

-  0.325 

0.326 

-  0.354 

0.038 

-  0.252 

0.035 

0.328 

-  0,415 

0.272 

-  0.379 

0.039 

-  0.397 

0.008 

0.220 

-  0.188 

0.257 

-  0.237 

0.048 

-  0.169 

0.048 

0.287 

-  0.269 

0.139 

-  0.326 

-  0.020 

-  0.115 

0.024 

0.141 

-  0.246 

0.051 

-  0.166 

-  0.089 

-  0.190 

-0.058 

0.099 

-  0.153 

0.123 

-  0.142 

-  0.004 

0.049 

0.114 

0.042 

-  0.169 

0.038 

-  0.103 

0.102 

-  0.129 

0.036 

0.102 

-  0.081 

0.113 

-  0.080 

0.028 

-  0.089 

0.012 

0.099 

-  0.116 

0.037 

-  0.126 

0*007 

-  0.030 

0.065 

0 

-73.407 

-10.096 


-  6.502 

-  3.468 

-  3.063 

-  1.399 

-  0.078 
0.186 
0.053 

-  0.388 

-  0.307 

-  0.246 

-  0.340 


Table  13 

Numerical  Values  of  the  Coefficients 
Aj[  and  Bi  (Equation  1)  for  Test  157 


Figure  1.  Typical  load  deformation  curves  obtained  from  the  direct 

shear  device  (site  7A) 


33 


Figure  2.  Peak  shear  load  versus  normal  load  (site  7A, 
test  series  98-102) 
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WET  DENSITY  =  100.4  PCF 
WATER  CONTENT  =64.5  PERCENT 


NORMAL  LOAD,  LB 

Peak  shear  load  versus  normal  load  (site  7B, 
test  series  107-111) 


Figure  4.  Peak  shear  load  versus  normal  load  (site  7A, 
test  series  117-122) 
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Figure  5.  Peak  shear  load  versus  normal  load  (site  7A, 
test  series  131-134) 
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Figure  6.  Peak  shear  load  versus  normal  load  (site  7B, 
test  series  139-144) 
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Figure  7.  Peak  shear  load  versus  normal  load  (site  8, 
test  series  157-159) 
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Figure  12.  Comparison  of  experimental  static 
failure  envelope  with  model  beha.'ior  (site  7B, 
test  series  i.39-144) 
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Figure  13.  Comparison  of  experimental  static 
failure  envelope  with  model  behavior  (site  8, 
test  series  157-159) 
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Figure  14.  Shear  stress-shear  displacement  relation;  field 
measurements  versus  model  behavior  (test  series  98-102); 
a  =  4.11  psi,  A  =  0.75  in. /sec 
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Figure  15.  Shear  stress-shear  displacement  relation;  field 
measurements  versus  model  behavior  (test  series  107-111) ; 
a  =  4.11  psi,  A  =  0.75  in. /sec 
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Figure  16.  Shear  stress-shear  displacement  relation;  field 
measurements  versus  model  behavior  (test  series  117-122); 
a  =  4.11  psi,  A  =  0.75  in. /sec 
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Figure  17.  Shear  stress-shear  displacement  relation;  field 
measurements  versus  model  behavior  (test  series  131-134); 
a  =  4.11  psi,  A  =  0.75  in. /sec 
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Figure  18.  Shear  stress-shear  displacement  relation;  field 
measurements  versus  model  behavior  (test  series  139— 14 A ) ; 
o  =  4.11  psi,  A  =  0.75  in. /sec 
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Figure  19.  Shear  stress-shear  displacement  relation;  field 
measurements  versus  model  behavior  (test  series  157-159); 
a  =  4.11  psi,  A  =  0.75  in. /sec 
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Figure  20.  Schematic  diagram  illustrating  the  direct  and 
indirect  field  measurements  for  a  circular  turn  test 


Figure  21.  Schematic  diagram  illustrating  the  validation 
procedure  for  the  WES  terrain-vehicle  interaction  model 
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Figure  22.  Outer  and  inner  track  velocity-time  histories  for 
test  99;  field  measurement  and  filtered  data  (input) 
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Figure  23.  Vehicle  speed-time  history  for  test  99;  comparison 
of  model  predictions  with  experimental  data 
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Figure  26.  Trajectory  of  the  center  of  gravity  of  the  vehicle 
for  test  99;  comparison  of  model  predictions  with  experimental 

data 
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Figure  27.  Outer  and  inner  track  velocity-time  histories  for 
test  100;  field  measurement  and  filtered  data  (input) 
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Figure  31.  Trajectory  of  the  center  of  gravity  of  the  vehicle 
for  test  100;  comparison  of  model  predictions  with  experimental 

data 
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Figure  32.  Outer  and  inner  track  velocity-time  histories  for 
test  101;  field  measurement  and  filtered  data  (input) 
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Figure  33.  Vehicle  speed-time  history  for  test  101;  comparison 
of  model  predictions  with  experimental  data 
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Figure  34.  Lateral  acceleration-time  history  for  test  101; 
comparison  of  model  predictions  with  experimental  data 
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Figure  37.  Outer  and  inner  track  velocity-time  histories  for 
test  102;  field  measurement  and  filtered  data  (input) 


TIME  T,SEC 


Figure  38.  Vehicle  speed-time  history  for  test  102;  comparison 
of  model  predictions  with  experimental  data 
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Figure  39.  Lateral  acceleration-time  history  for  test  102; 
comparison  of  model  predictions  with  experimental  data 
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Figure  41. 
test  102; 


Trajectory  of  the  center  of  gravity  of  the  vehicle  for 
comparison  of  model  predictions  with  experimental  data 
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Figure  42.  Outer  and  inner  track,  velocity-time  histories  for 
test  107;  field  measurement  and  filtered  data  (input) 
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Figure  44.  Lateral  acceleration-time  history  for  test  107; 
comparison  of  model  predictions  with  experimental  data 
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Figure  45.  Total  power-time  history  for  test  107;  comparison 
of  model  predictions  with  experimental  data 
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Figure  46.  Trajectory  of  the  center  of  gravity  of  the  vehicle  for 
test  107;  comparison  of  model  predictions  with  experimental  data 
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Figure  48.  Vehicle  speed-time  history  for  test  108;  comparison 
of  model  predictions  with  experimental  data 
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Figure  52.  Outer  and  inner  track  velocity-time  histories  for 
test  109;  field  measurement  and  filtered  data  (input) 
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Figure  55.  Tocal  power-time  history  for  test  109;  comparison 
of  model  predictions  with  experimental  data 
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Figure  57.  Outer  and  inner  track  velocity-time  histories  for 
test  110;  field  measurement  and  filtered  data  (input) 


Figure  53.  Vehicle  speed-time  history  for  test  110;  comparison 
of  model  predictions  with  experimental  data 
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Figure  62.  Outer  and  inner  track  ve.Locity-time  histories  for 
test  111;  field  measurement  and  filtered  data  (input) 
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Figure  64.  Lateral  acceleration-time  history  for  test  111; 
comparison  of  model  predictions  with  experimental  data 
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Figure  65.  Total  power-time  history  for  test  111;  comparison 
of  model  predictions  with  experimental  data 
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Figure  67.  Outer  and  inner  track  velocity-time  histories  for 
test  117;  field  measurement  and  filtered  data  (input) 
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Figure  69.  Lateral  acceleration-time  history  for  test  117; 
comparison  of  model  predictions  with  experimental  data 
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Figure  71.  Trajectory  of  the  center  of  gravity  of  the  vehicle 
for  test  117;  comparison  of  model  predictions  with  experimental 

data 
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Figure  72.  Outer  and  inner  track  velocity-time  histories  for 
test  118;  field  measurement  and  filtered  data  (input) 
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Figure  73.  Vehicle  speed-time  history  for  test  118;  comparison 
of  model  predictions  with  experimental  data 
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Figure  77.  Outer  and  inner  track  velocity-time  histories  for 
test  133;  field  measurement  and  filtered  data  (inp  ) 
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Figure  78.  Vehicle  speed-time  history  for  test  133;  comparison 
of  model  predictions  with  experimental  data 
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Figure  80.  Total  power-time  history  for  test  133;  comparison 
of  model  predictions  with  experimental  data 
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Figure  82.  Outer  and  inner  track  velocity-time  histories  for 
test  134;  field  measurement  and  filtered  data  (input) 
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Figure  83.  Vehicle  speed-time  history  for  test  134;  comparison 
of  model  predictions  with  experimental  data 
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Figure  85,  Total  power-time  history  for  test  134;  comparison 
of  model  predictions  with  experimental  data 
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Figure  86.  Trajectory  of  the  center  of  gravity  of  the  vehicle  for 
test  134;  comparison  of  model  predictions  with  experimental  data 
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Figure  88.  Vehicle  speed-time  history  for  test  139;  comparison 
of  model  predictions  with  experimental  data 


LATERAL  ACCELERATION 


115 


525 


LEGEND 


-  FIELD  MEASUREMENTS 


0  20  40  60  80  1 00 

TIME  T,SEC 


Figure  90.  Total  power-time  history  for  test  139;  comparison 
of  model  predictions  with  experimental  data 
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gure  91.  Trajectory  of  the  center  of  gravity  of  the  vehicle  fi 
;est  139;  comparison  of  model  predictions  with  experimental  dat. 
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Figure  93.  Vehicle  speed-time  history  for  test  140;  comparison 
of  model  predictions  with  experimental  data 
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Figure  95.  Total  power-time  history  for  test  140;  comparison 
of  model  predictions  with  experimental  data 
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Figure  102.  Outer  and  inner  ““^jJ^ered^at^Unplt) 
test  157;  field  measurement  and  filtere 
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Figure  103.  Vehicle  speed-time  history  for  test  157;  comparison 
of  model  predictions  with  experimental  data  (example  of  poor 

quality  data) 
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Figure  105.  Total  power-time  history  for  test  157; 
comparison  of  model  predictions  with  experimental 
data  (example  of  poor  quality  data) 
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Figure  106.  Trajectory  of  the  center  of  gravity  of  the  vehicle  for 
test  157;  comparison  of  model  predictions  with  experimental  data 
(example  of  poor  quality  data) 


Figure  107.  Turning  radius  versus  steering  ratio  for  test  series 
107-111;  comparison  of  model  predictions  with  experimental  data 
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Figure  108,  Inner  track  velocity  versus  turning  radius  for  test  series  107-111;  comparison  of 
model  predictions  with  experimental  data  (for  tests  at  radii  of  33.74  and  63.69  ft,  stability 

criteria  govern  the  speed  of  the  vehicle) 
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Figure  109.  Outer  track  velocity  versus  turning  radius  for  test  series  107-111; 
comparison  of  model  predictions  with  experimental  data 
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Figure  110.  Vehicle  speed  versus  turning  radius  for  test  series  107-111; 
comparison  of  model  predictions  with  experimental  data 
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Figure  111.  Power  requirement  versus  turning  radius  for  test  series  107-111; 
comparison  of  model  predictions  with  experimental  data 
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Figure  112.  Lateral  acceleration  versus  turning  radius  for  test  series  107-111; 
comparison  of  model  predictions  with  experimental  data 


APPENDIX  A:  TERRAIN-VEHICLE  INTERACTION  MODEL 


Introduction 


1.  The  basic  concepts  of  the  theory  of  terrain-vehicle  interaction 
were  developed  during  the  1950's  by  Bekker  (1963).*  By  assuming  various 
load  distributions  along  the  tracks,  Bekker  was  able  to  develop  several 
mathematical  expressions  relating  the  characteristics  of  the  vehicle  and 
the  tractive  effort  of  the  terrain  during  steering,  ly  considering  the 
lateral  and  longitudinal  coefficients  of  friction  between  the  track  and 
the  ground,  Hayashi  (1975)  developed  simple  equations  for  practical 
analysis  of  steering  of  tracked  vehicles.  Hayashi's  work,  however,  did 
not  include  the  effect  of  the  centrifugal  forces  on  steering  performance 
of  the  vehicle.  Kitano  and  Jyorzaki  (1976)  developed  a  more  compre¬ 
hensive  model  for  uniform  turning  motion  including  the  effects  of 
centrifugal  forces.  This  model,  however,  is  based  on  the  assumption 
that  ground  pressure  is  concentrated  under  each  road  wheel  and  the 
terrain-track  interaction  is  simulated  by  Coulomb-type  friction.  The 
model  given  by  Kitano  and  Jyorzaki  was  extended  by  Kitano  and  Kuma 
(1977)  to  include  nonuniform  (transient)  motion,  but  the  basic  elements 
of  the  terrain-track  interaction  part  of  the  model  were  retained. 

Baladi  and  Rohani  (1978)  developed  a  model  for  uniform  turning  motion 
parallel  to  the  development  by  Kitano  and  Jyorzaki  insofar  as  the 
kinematics  of  the  vehicle  are  concerned.  In  contrast  to  the  development 
by  Kitano  and  Jyorzaki  (1976),  however,  this  model  is  based  on  a  more 
comprehensive  soil  model.  Baladi  and  Rohani  (1979  and  1981)  extended 
the  WES  terrain-vehicle  model  completed  in  1979  to  include  nonuniform 
(transient)  motion  on  level  terrain.  In  addition,  the  WES  soil  model 
was  modified  to  include  a  nonlinear  failure  envelope  describing  the 
shearing  strength  of  the  terrain  material  (Baladi  and  Rohani,  1981). 


*  References  cited  in  this  appendix  are  listed  in  the  References  sec¬ 
tion  at  the  end  of  the  main  text. 
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In  this  appendix,  the  WES  terrain-vehicle  model  is  extended  to  include 
the  treatment  of  nonuniform  (transient)  turning  motion  on  sloping 
terrain. 


Soil  Model 


Strength  components 

2.  One  of  the  most  important  properties  of  soil  affecting  traffic- 
ability  is  in  situ  shear  strength.  It  has  been  found  experimentally 
that  the  shear  strength  of  purely  cohesive  soils  is  relatively  indepen¬ 
dent  of  the  confining  stress,  but  strongly  affected  by  the  time  rate  of 
shearing.  On  the  other  hand,  the  shear  strength  of  purely  frictional 
soils  is  found  to  be  relatively  independent  of  the  time  rate  of  loading, 
but  strongly  dependent  on  confining  pressure.  The  shearing  resistance 
of  most  soils,  however,  is  due  to  both  the  frictional  and  cohesive 
components.  The  cohesive  and  frictional  components  of  strength  are 
usually  added  together  in  order  to  obtain  the  total  shear  strength  of 
the  material.  For  static  loading  (very  slow  rate  of  deformation),  the 
shear  failure  envelope  is  defined  by 

=  A  -  M  exp(-Na)  (Al) 

where 

=  the  maximum  shearing  strength  of  the  material 
a  =  the  normal  stress 

A  =  the  strength  of  the  material  when  o  is  large 
A-M  =  C  =  the  strength  of  the  material  or  cohesion  when  a  =  0 
N  =  a  material  constant 

Equation  Al  is  shown  graphically  in  Figure  Al. 

3.  As  noted  previously,  the  shear  strength  of  cohesive  soils 
increases  with  the  increasing  rate  of  loading.  For  the  range  of  load¬ 
ing  rates  associated  with  the  motion  of  tracked  vehicles,  the  contribu¬ 
tion  to  cohesive  strength  due  to  dynamic  loading  can  be  expressed  as 
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C,[l  -  exp(-AA)]  ,  where  C,  and  A  are  material  constants  and  A  is 
d  d 

the  time  rate  of  shearing  deformation.  In  view  of  the  above  expres¬ 
sion,  the  dynamic  failure  criterion  takes  the  following  form: 

=  A  +  C^[l  -  exp(-AA)]  -  M  exp(-No)  (A2) 


When  A  equals  zero,  the  dynamic  failure  criterion  (Equation  A2) 
reduces  to  the  static  failure  criterion  (Equation  Al).  Both  are  shown 
graphically  in  Figure  Al. 

Shear  stress-shear  deformation  relation 

4.  Prior  to  failure,  the  shear  stress-shear  deformation  charac¬ 
teristics  of  a  variety  of  soils  can  be  expressed  by  the  following 
mathematical  expression  (Kondner,  1963): 


Gtma 
tm  +  GIAI 


(A3) 


The  behavior  of  Equation  A3  is  shown  graphically  in  Figure  A2,  in  which 
T  denotes  shearing  stress,  A  is  shearing  deformation,  and  G  is  the 
initial  shear  stiffness  coefficient.  Substituting  xM  from  Equation 
A2,  the  shear  stress-shear  deformation  relation  for  soil  becomes 


G[A  +  C ,  -  C ,  exp(-AA)  -  M  exp (-No)] A 

T  =  - - - - - : -  (A4) 

G |  A  |  +  A  +  exp(-AA)  -  M  exp (-No) 

For  purely  cohesive  soils,  N  equals  zero  and  t  is  only  a  function  of 
A  and  A  .  For  cohesionless  or  granular  soils,  M  equals  A  , 

is  zero,  and  t  is  a  function  of  A  and  a  .  For  mixed  soils 
exhibiting  shearing  resistance  due  to  both  frictional  and  cohesive 
components,  x  is  dependent  on  A  ,  A  ,  and  a  .  The  qualitative 
behavior  of  Equation  A4  for  these  three  conditions  is  shown  in  Figure 
A3.  It  should  be  pointed  out  that  Equation  A4  reduces  to  the  rigid 
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plastic  soil  model  often  used  in  mobility  studies  when  an  extremely 
large  value  is  specified  for  G  and  A  and  is  set  to  zero. 

5.  An  appropriate  test  for  determining  the  numerical  values  of 
the  six  material  constants  in  Equation  A4  is  an  in  situ  direct  shear 
test.  A  field  direct  shear  device  has  been  developed  at  the  WES  for 
this  purpose.  A  description  of  this  device  and  the  method  of  analysis 
of  the  data  obtained  from  the  direct  shear  test  are  documented  in 
Appendix  C  of  a  report  by  Baladi  and  Rohani  in  1979.  For  completeness, 
this  Appendix  C  is  included  herein  as  Appendix  B. 

Derivation  of  Terrain-Vehicle  Model 


Boundary  condition 

6.  The  geometry  of  the  vehicle  and  the  boundary  conditions  of  the 
proposed  model  are  shown  schematically  in  Figure  A4.  The  XYZ  coordi¬ 
nates  are  the  local  coordinate  system  of  which  X  is  always  the  longi¬ 
tudinal  axis  of  the  vehicle  and  Y  is  a  transverse  axis  parallel  to 
the  ground.  These  axes  intersect  at  the  center  of  geometry  of  the 
vehicle  0  .  The  Z  axis  is  a  vertical  axis  passing  through  the 
origin  0  .  The  center  of  gravity  of  the  vehicle  (CG)  lies  on  the 

X  axis  and  is  displaced  by  a  distance  from  the  origin.  The 

numerical  value  of  is  assumed  to  be  positive  if  CG  is  displaced 

forward  from  the  center  of  geometry  of  the  vehicle.  The  XY  coordi¬ 
nates  of  the  instantaneous  center  of  rotation  (ICR)  are  P  +  C^  and 
R  ,  respectively,  where  P  is  the  offset.  The  center  of  rotation  and 
the  radius  of  the  trajectory  of  the  CG  are,  respectively,  CR  and 
R0  .  The  height  of  the  center  of  gravity  measured  from  ground  surface 
is  denoted  by  H  .  The  lengths  of  the  track-ground  contact,  the  track 
width,  and  the  tread  of  the  tracks  are  L  ,  D  ,  and  B  ,  respectively. 
As  shown  in  Figure  A4,  the  components  of  the  inertial  forces  Fq 
in  X  and  Y  directions  are,  respectively,  F^  and  F^y  .  The 
weight  of  the  vehicle  is  W  . 

Stress  distribution  along  the  tracks 

7.  Two  types  of  stress  (i.e.,  normal  and  shear  stresses)  exist 
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along  the  track.  As  indicated  in  Figure  A4,  th*  normal  stresses  under 
the  outer  and  inner  tracks  are  denoted  by  R^(X)  and  [^(X)  ,  respec¬ 
tively.  The  components  of  the  shear  stress  in  X  and  Y  directions 
are,  respectively,  T^(X)  and  Q^(X)  for  the  outer  track,  and  T2(X) 
and  (^(x)  f°r  the  inner  track.  These  stresses  are  dependent  on  the 
terrain  type,  vehicle  configuration,  and  speed  and  turning  radius  of 
the  vehicle. 

8.  The  magnitude  of  normal  stresses  R^(X)  and  R2W  can  be 
determined  in  terms  of  the  components  of  the  inertial  force,  the  track 
tensions,  and  the  characteristics  of  the  vehicle  by  considering  the 
balance  of  vertical  stresses  and  their  moments  in  Figure  A4.  Thus 


R1  (x)  =  ~^r 
1  dl/  L 


1  .  ,  h  fcy  Fcx  .  dL  N1(X) 

2  +  6xcx  -  b  IT  -  6hx  IT  +  — w~~ 


(A5) 


R?(x)  =  -~ 
dL 


1  h  FCY  FCY  dL  N2(x) 

2  +  6XCX  +  b  ~  '  6hx  IT  +  — H - 


(A6) 


where 

h  =  H/L 
b  =  B/L 
d  =  D/L 

cx  =  Cx/L 

x  =  X/L 
y  =  Y/L 
z  =  Z/L 

N^(x)  and  ^(x)  =  contributions  due  to  track  tension. 


9.  The  components  of  the  shear  stress  in  the  X  and  Y  direc¬ 
tions  along  both  the  outer  and  inner  tracks  can  be  obtained  by  combining 
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Equations  A4,  A5,  and  A6.*  Thus  (it  is  noted  that  and  replace 

the  normal  stress  a  in  Equation  A4) 


T  (x)  = 

L 


da  +  dCj  -  dc^  exp(-A<5^)-m  exp[-nr^(x)] 


6^|d  +  da  +  dc^  -  dc^  exp (-A6^) -m  exp[-nr^(x)] 


cosy^  (A7) 


Q  (x)  =  ^yy5, 
L 


da  +  dc^  -  dc^  exp(-A6^)-m  exp[-nr^(x)] 


6.  d  +  da  +  dc,  -  dc,  exp(-A6.)-m  exp[-nr.(x)] 
1  d  d  l  r  i 


sinY.  (A8) 
i 


where 


i  =  1,2 


r^(x)  =  dL  R^(x)/W 

6.  =  A./L 
1  x 


si  ’  VL 


y  =  GL  /W 
A  =  AL 

a  =  AL2/W 

m  =  ML2/W 

n  =  NW/L2 


Cd  =  CdL  /W 

The  variables  Y^  and  Y ^  >  in  Equations  A7  and  A8,  are  the  slip 
angles  and  can  be  written  as 


*  To  account  for  the  effect  of  the  size  of  the  shear  box  on  the  shear 
stiffness  G  ,  the  measured  value  of  G  is  normalized  by  multiplying 
it  by  4/L  (the  length  of  the  shear  box  =  4  in.). 
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where 


P  =  P/L 

The  parameter  is  the  distance  between  the  instantaneous  center  of 

rotation  of  the  outer  track  IC^  and  its  axis  of  symmetry,  and  C2  is 
the  distance  between  the  instantaneous  center  of  rotation  of  the  inner 
track  IC2  and  its  axis  of  symmetry  (Figure  A5) . 

10.  In  order  to  use  Equations  A7  through  A9,  the  normal  stress 

contributions  due  to  track  tensions  N^(x)  and  ^(x)  ,  the  track  slip 

•  • 

velocities  and  displacements  (i.e.,  ,  A^  ,  ,  and  ,  and  the 

inertial  forces  F  and  F  have  to  be  determined.  These  factors 
LX  Lx 

are  discussed  in  the  following  paragraphs. 

Normal  stress  contribution  due  to  track  tension 

11.  The  effect  of  track  tension  on  the  normal  stress  distribution 
is  influenced  considerably  by  the  motion  of  the  vehicle.  At  relatively 
low  speed,  tractive  effort  is  applied  to  the  outer  track,  while  braking 
force  is  applied  to  the  inner  track  (Figure  A6a) .  At  high  speed,  on 
the  other  hand,  tractive  efforts  are  applied  to  both  tracks  (Figure 
A6b) . 

12.  The  angles  0  and  9,  in  Figure  A6  are  the  approach  and 
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departure  angles  of  the  track  envelope,  respectively.  The  forces  T^ 
and  T2  are  the  track  tension  in  the  outer  and  inner  track,  respective¬ 
ly.  These  forces  can  be  obtained  by  integrating  Equation  A7.  Thus 


A7 


T 


1 


1 

2 


f' 


(x)  dx 


T 


2 


(A10) 


The  normal  stress  distributions  are  influenced,  however,  by  the  vertical 
components  of  the  forces  and  T2  ;  namely,  n^  ,  n^  ,  and  n^  . 

The  values  of  ,  n2  ,  and  are 


n,  =  T.,  sin  0 , 
11  d 


n2  ' 


T„  sin  0  if  ?0  >  0 

4  H  L 


if  S2  <  0 


(All) 
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-T,,  sin  0d  if  ?2  <  0 


if  ?2>o 


(A13) 


With  the  determination  of  the  forces  ,  n2  ,  and  n2  ,  the  normal 
stress  contributions  due  to  track  tension  may  be  determined. 

13.  Since  the  tracks  are  assumed  to  be  rigid,  the  normal  stresses 
due  to  track  tension  may  be  distributed  according  to  the  following 
equations  (Figure  A7): 
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N,  (x)  =  ax  +  m 
1  o 


,  *  1  ,  „  1 
for  l"2^X^2 


2n 


„  ,  v  .  .  1  /  .  1  t\  -  1  „  i  1 

N,  (x)  =  ax  +  m  +  — r  x  +  -r  -  —  for  -  -r  <  x  <  - —  -= 

1  0  cU2  1  2  L/  2  -  -L  2 


(A14) 


and 


„  ,  ,  .  2n2  /  1  .  I \  ,  111 

N2(x)  =  ax  +  m1  -  -2  x  -  y  +  J  f°r  2_L-X-2 

dx. 


N2(x)  =  ax  + 


,  *  1  „  „  1  l 

for  r2^XirL 


2ni 


2  /  l  £\  1  £  1 

N2(x)  ~  ax  +  irij  +  — ^  ^x  +  -j  -  -j  for  -  y  £  x  ^  —  -  -j 

di 


(A15) 


in  which  i  is  the  distance  between  two  adjacent  wheels,  and  a  , 
idq  ,  and  can  be  determined  by  considering  the  equation  of  equilib¬ 

rium  of  normal  stresses  and  the  moments  of  these  stresses.  Thus 


i  1 


2  L  f  2  2n 

(  (ax  +  m  )dx  +  I  — ir  (x  +  \  -  7 )  dx  =  0 

Jl  °  *1  dZ2  1  2  L/ 


(A16) 


dx 


2  l, 

I  1 

L  “  2 


[  2^1  /  +1  A\ 

J  d£2  (X  2_L^ 


dx  =  0 


(A17) 
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and 


(2ax  +  m 


+  m  ) 
o 


dx 


dx 


(A18) 


1 

f  2(n1  +  n'2) 

+  Jl  d£2 
2 


dx  =  0 


Equations  A16  through  A18  contain  three  unknowns:  a  ,  mQ  >  an^  mj  • 
Completing  the  integrations  results  in 


a 


-  n' 


(A19) 


m 

o 


dL 


2  "l 


(A20) 


=  —^2  ^  +  np  (A21) 

dL 


Substitution  of  Equations  A19  through  A21  into  Equations  A14  and  A15 
leads  to 

Nx(x) 

N^(x) 

for  -  y  £  x  £  6  -  y 


=  —^2  [(3  -  26)  (n2  “  n2  “  n^x  +  for  6  -  y  £  x  £  y 


dL 

1 

dL2 


2n 

(3  -  26)  (n  -  n'  -  n .)  +  ~y 
*  6 


x  +  (i 


-  6 


(A22) 
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and 


N„(x)  =  — ^ 
dL 


(3  -  28)  (n2  ~  n2  "  n^)  - 2~ 


rt 

x  +  n2  +  n' 


for  -  8  <  x  <  j 


N2(x)  =  — 2  [(3  -  28)  (n2  -  n2  -  x  +  n2  +  n2] 
dL 


(A23) 


for  8-'|'£xf_y-8 


N„(x)  =  — j 
dL 


2nl 

(3  -  28)  (n2  -  -  nx)  +  — 


X  +  (H-8)  "2  +  "2 


for  ~  x  <_  8  - 


where 


8  =  £  (A24) 

Note  that  Equations  A12  and  A13  dictate  that  either  n2  or  n2  in 
Equations  A22  and  A23  is  zero. 

Kinematics  of  the  vehicle 

14.  A  tracked  vehicle  in  transient  motion  is  shown  schematically 

in  Figure  A8.  The  XYZ  coordinates  are  the  local  coordinate  systems 

that  are  fixed  with  respect  to  the  moving  vehicle  (also  see  Figure  A4). 

The  origin  0  of  this  coordinate  system  stays,  for  all  time,  at  a 

distance  C  from  the  center  of  gravity  of  the  vehicle.  The  ¥$ 

A 

coordinate  system  is  fixed  on  level  ground,  and  its  origin  coincides 


All 


with  the  center  of  gravity  at  time  zero.  The  vehicle  can  maneuver  on 
the  Ti*  plane  and  the  displacements  of  the  center  of  gravity  of  the 
vehicle  from  this  reference  frame  are  T(t)  and  $(t)  . 

15.  The  velocities  v  and  v  (relative  to  the  origin  of  the 
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coordinate  system)  as  well  as  the  velocities  v^  and  v^  are 
related  to  the  instantaneous  velocity  v  of  the  CG  by 


=  VVX  +  VY  =  VVf  + 
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The  side-slip  angle  a  ,  which  is  the  angle  between  the  velocity  vector 

v  and  the  longitudinal  X  axis  of  the  vehicle,  is  related  to  the 

velocities  v  and  v  as 
A  JL 


-1  VY  da  /  dVY  dVx\  /  2 
a  =  tan  -  ,  ^  ^  -  vy  —jjv 


(A26) 


The  yaw  angle  u  and  the  directional  angle  0  are  related  to  a  as 


0  =  a)  -  a  , 


dO  _  da)  da 
dt  dt  dt 


(A27) 


Substitution  of  Equation  A26  into  Equation  A27  leads  to 


\ 


d9  dm 
dt  dt 


(A28) 


16.  The  radius  of  curvature  of  the  trajectory  of  the  center  of 
gravity  (i.e.,  the  distance  between  CR  and  CG  (Figures  A5  and  A9))  is 


R 

o 


v 


/d0 
/  dt 


dui 

dt 


3 

v 


v 


dvy  dv 

- 4*  v  — 

X  dt  Y  dt 


X 
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The  coordinates  of  the  trajectory  of  the  center  of  gravity  of  the 
vehicle  can  be  written  as 


(A30) 


17.  The  coordinates  of  the  instantaneous  center  of  rotation  (ICR) 
of  the  hull  in  the  XY  systems  (X^.  ,  Y^.)  and  the  instantaneous  radius 
of  curvature  (R^)  are  (Figures  A5  and  A9) 

XI-P  +  CX  =  ''Y /ff+CX 

ft  (A31> 

Rj  =  \]rZ  +  P2 


The  instantaneous  velocities  of  an  arbitrary  point  e  of  the  hull  are 
shown  in  Figure  A9  and  can  be  written  as 


eX 


=  v„ 


v  dW 
dt 


eY 


-  (X  - 


v£ 


(A32) 


v  = 
e 


do) 

dt 


2 
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Track  slip  velocity  and  displacement 


18.  Assume  that  (vg^  =  A^)  is  the  slip  velocity  of  an 

arbitrary  point  e^  of  the  outer  track  and  v^  ( vg ^  is  the 

slip  velocity  at  point  (e^  and  e2  have  the  same  abscissa)  of  the 
inner  track  (Figure  A5).  The  X  and  Y  components  of  these  velocities 
are 


_  dio  _  T  dm 
VsXl  "  C1  dt  h  L  dt 


Vl  “  <X  -  P  -  CX>  ft  -  L«  -  CX>  %  -  VY 


For  the  outer  track  (A33) 


v 


sX2 


*2 


dm 

dt 


v 


sY2 


sYl 


’  For  the  inner  track 

. 


(A34) 


As  indicated  in  Figure  A10,  the  angular  velocity  du)/dt  and  the  value 
of  R  can  be  written  as 


dm  _  1  . 

dt  =  b£  UX1  ”  VsXl 


X2 


+ 


\x2> 


R  "  tdS  (VX1 
dt 


VsXl  + 


X2 


-  vsX2) 
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where 


v  =  the  velocity  of  the  outer  track  in  X  direction 
XI 

VX2  =  t^e  ve^-oc:'-tY  the  inner  track  in  X  direction 

The  ratio  of  v^  and  v^  is  defined  as  the  steering  raHo  e  . 
Thus 


e 


VX1'/VX2 


(A36) 
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Substitution  of  Equations  A31  and  A36  into  Equation  A35  leads  to 


_  /  bL  du)\ 

VsXl  “  X2  (VX  +  2  dt) 

/  bL  dui\ 

VsX2  =  VX2  "  (VX  2  dt) 


For  the  outer  track 


For  the  inner  track 
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Comparison  between  Equations  A37  and  A38  and  Equations  A33  and  A34 
results  in 


(A38) 


The  slip  velocities  and  displacements  of  the  outer  and  inner  tracks  can 
be  obtained  from  Equations  A33,  A34,  and  A37.  Thus 


si 


dt  + 


II 

L 


(A39) 


(A40) 


\ 


where 

t:  =  (L/2  -  X)/vn 
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t2  =  (L/2  -  X)/vx2 


/ 


A^  =  initial  displacement  of  the  outer  track 

Aj2  =  initial  displacement  of  the  inner  track 

The  balance  of  forces  and  moments  dictates  that  these  initial  displace¬ 
ments  be  numerically  equal  to  is  the  coefficient  of  rolling 

resistance,  which  must  be  measured  experimentally  or  calculated  from 
empirical  relations  presented  by  Rula  and  Nuttall  (1971)  for  each  soil 
type  and  each  vehicle) . 

Inertial  forces 

19.  According  to  Figure  A8,  the  relationship  between  the  veloci¬ 
ties  Vy  and  v^  and  the  velocities  vx  and  v^  can  be  written  as 


v  =  -v  cos  a)  -  v  sin  ui 

i  X  I 

v.  =  v„  sin  !»)  -  v„  cos  m 
$  X  Y 

j 

The  acceleration  in  T  and  <t>  direction,  a^  and  a$  ,  can  be 
written  as 


(A41) 


aT  "  dt 


a$  dt 


(A42) 


The  forward  and  lateral  accelerations,  ax  and  a^  ,  can  be  written  in 
terms  of  a^  and  a^  as 


a,,  =  -a,,,  cos  a)  +  a.  sin  to 
X  V  $ 

a^  =  -a^  sin  u>  -  a^  cos  w 


(A43) 


\ 


t 
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Substitution  of  Equations  A41  and  A42  into  Equation  A43  leads  to 


and 


_  du 

3X  "  dt  +  VY  dt 


_  S  du 

S  dt  VX  dt 


(A44) 


Hence,  the  X  and  Y  components  of  the  inertial  force  can  be  written 
as 


v  _w  _W[^X  du] 

Fcx  g  aX  g  l dt  VY  dt/ 


and 


W  H  /  dVY 


du 


FCY  ~  g  ®Y  “  g  l  dt  "  VX  dt 


(A45) 


The  rolling  resistance 

20.  The  rolling  resistance  Rs  is  a  function  of  terrain  type, 
vehicle  speed,  track  condition,  etc.  Therefore,  rolling  resistance 
should  be  measured  for  every  specific  condition.  In  this  formulation, 
however,  the  rolling  resistance  is  assumed  to  be  proportional  to  normal 
load .  Thus 

1 

} 

R  =  — r  /  [r.  (x)  +  r?(x)]dx  (A46) 

S  dL  A  1 

"  2 


Equation  of  motion 

21.  Steerability  and  stability  of  tracked  vehicles  depend  on  the 
dynamic  balance  between  all  forces  and  moments  applied  on  the  vehicle. 
According  to  Figure  A4,  the  following  three  equations  govern  the  motion 
of  the  vehicle: 
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(A47) 


j  [t1(x)  +  t2(x)]  dx  -  i  [^(x)  +  r2(x)]dx  -  f 


CX 


j  [qL(x)  +  q2(x)]  dx  =  f 


CY 


(A48) 


(x)  +  q,(x)]  (x  -  c  )  dx  + 


Utei 


(x)  - 


t2(x)]  dx 


I 


^  i.  d  **03 

+  |  i  J  “  rl(x)^  dx  =  LW  ^2 


(A49) 


where 


•!«  '  T1W 
‘  ¥  T2W 
-  t-  Qi« 

dL2 


q2(x)  =  IT  Q2(x) 


:  =  _CX 

Ccx  w 

E 

CCY  W 


morrent  of  inertia  about  an  axis  passing  through  the 
^  and  oarallel  to  the  Z  axis  (Figure  A4) 


X  — •  mass  _ _ _ 

z  CG  of" the  vehicle  and  parallel 
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Equations  A47  through  A49  with  the  aid  of  Equation^  A7  through  A46 
constitute  three  equations  that  involve  three  unknow.  s.  The  three 
unknowns  are  either  v^  ,  ,  and  drn/dt  or  ^  »  anc*  P  • 

order  to  obtain  a  complete  solution  for  either  of  the  two  sets  of 
unknowns,  one  of  the  following  driving  conditions  must  be  specified: 

(a)  time  history  of  the  steering  ratio  e ( t )  and  the  initial  speed  of 
the  vehicle,  (b)  time  history  of  the  velocity  of  the  individual  tracks 
v^(t)  and  v^(t)  an^  t*ie  initial  speed  of  the  vehicle,  (c)  time 
history  of  the  velocity  of  the  vehicle  v(t)  and  the  trajectory  of 
motion,  (d)  time  history  of  the  velocity  of  the  vehicle  and  a  constant 
value  of  steering  ratio  e  ,  or  (e)  the  trajectory  of  motion  and  a 
determination  of  the  maximum  velocity  time  history  at  which  the  vehicle 
can  traverse  the  specified  trajectory.  A  computer  program  called  AGIL 
was  developed  to  solve  Equations  A47  through  A49  using  Newton's  itera¬ 
tion  technique. 

Treatment  of  sloping  terrain 

22.  The  analysis  of  the  nonuniform  turning  motion  on  sloping 
terrain  presented  in  this  section  is  an  extension  of  the  corresponding 
analysis  for  steady-state  motion  on  sloping  terrain  reported  by  Baladi 
and  Rohani  (1979  and  1981).  Figure  All  shows  schematically  a  tracked 
vehicle  under  nonuniform  (transient)  turning  motion  on  a  terrain  with 
slope  angle  q  .  In  this  case,  the  weight  of  the  vehicle  W  could  be 
resolved  into  a  normal  component  (normal  to  the  terrain)  and  a 

parallel  component  .  Thus 


=  W  cos  q 

- 

WT  =  W  sin  n. 


(A50) 


In  general,  the  longitudinal  axis  of  the  vehicle  X  makes  an  angle  x 
with  the  component  (Figure  All).  Therefore,  the  component 

could  be  resolved  into  two  components.  The  first  component  W  is 

1a 


parallel  to  the  X  axis  of  the  vehicle  and  the  second  component 
is  parallel  to  the  Y  axis.  Thus 


WTX  =  WT  cos  X  =  W  sin  n  cos  x 


(A51) 


WTY  =  WT  sin  x  =  W  sin  n  sin  x  J 

The  angle  x  is  related  to  the  yaw  angle  w  through  the  following 
relation 


X  =  m  +  v  (A52) 

where  v  is  a  constant.  The  numerical  value  of  v  depends  on  the 
initial  position  of  the  vehicle  (v  =0,  90,  180,  and  270  deg  corre¬ 
sponds,  respectively,  to  the  initial  position  of  the  vehicle  at  points 
1,  2,  3,  and  4  in  Figure  All).  Substitution  of  Equation  A52  into 
Equation  A51  leads  to 


W 


TX 


=  W  sin  n  cos (w  +  v) 


WTY  =  W  sin  n  sin(w  +  v)  J 


(A53) 


In  view  of  Equations  A50  and  A53,  the  normal  stresses  under  the  outer 
and  inner  tracks  (Equations  A5  and  A6)  become 


R,  (x)  = 

dL 


cos  n  .  ,  h 

2  +  6xcx  cos  n  -  b 

CY 

— - sin  n  sin(a)  +  v) 

-  6hx 


CX 

— -  +  sin  q  cos(o)  +  v) 


DLNx(x) 


W 
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n  /  \  w  J cos  n  .  ,  .  h 

R2(x)  =  — 2  ' — 2 —  +  6xCX  C0S  n  +  b’ 
dL  l 


CY 

L  w 


-  sin  n  sin(w  +  v) 


6hx 


CX 


+  sin  n  cos(w  +  v) 


DLN2(x)' 
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Equations  A54  and  A55  can  be  combined  with  Equations  A7  through  A46  to 
develop  the  equations  of  motion  for  a  sloping  terrain.  Thus 


i  i. 

/[tl(x)  +  t2(x)]  dx  -  J  [rx(x)  +  r2(x)]  dx 


=  fcx  +  sin  n  cos  (id+v) 


(A56) 


f  [q-j^x)  +  q2(x)]  dx  =  fCY  -  sin  n  sin(u  +  v) 


(A57) 


f  U^x)  +  q2^x^  ^x  "  cx)  dx  +  \  f  [t2(x)  -  t^x)]  dx 


+  2  6  f  [ri(x)  -  r2(x)]  dx  =  lt7T 

dt 


(A58) 


A21 


Sprocket  power 

23.  The  steering  performance  of  a  tracked  vehicle  may  be  limited 
either  by  its  stability  or  by  the  power  available  at  the  sprockets. 

The  powers  that  must  be  available  at  the  inner  and  outer  track 
sprockets,  PT1  and  PT2  ,  respectively,  are 


~  f  [t,  (x)  6  /cos  yj 

w/EF  /lF  4  1  1  1 


dx 


I  r  (x)  dx 
J1 


XI  2 
+  a  cos  n - ; - 


f  [rx(x)  +  r2(x)J  dx 


(A59) 


PT2 

W/Lg 


~  /  tc2(x) 

/Eg  J1 

2 


S^/cos  y2]  dx 


I  r„(x)  dx 
•'l 


i  X2  2 

+  |5  cos  n - r 


i 


f  [r1(x)  +  r2(x)]  dx 


(A60) 


A22 


Therefore,  the  total  power 
required  are 


PT  and  the  differential  power  PTD 


PT  =  PT1  +  PT2 
PTD  =  PT1  -  PT2J 


(A61) 
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Figure  Al.  Proposed  failure  relation  for 


SHEARING  STRESS, 


Figure  A2.  Proposed  shear  stress/deformation  relation 
during  shearing  process  (Equation  A3) 
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SHEARING  STRESS,  t  SHEARING  STRESS,  r  SHEARING  STRESS, 


SHEARING  DEFORMATION,  A 

a.  PURELY  COHESIVE  SOIL  (r  INDEPENDENT  OF  <r> 


SHEARING  DEFORMATION,  A 

b.  GRANULAR  SOIL  (r  INDEPENDENT  OF  A) 


SHEARING  DEFORMATION,  A 

c.  MIXED  SOIL  (t  DEPENDENT  ON  BOTH  a  AND  A) 


Figure  A3.  Qualitative  behavior  of  the  soil 
model  (Equation  A)  for  various  types  of  soil 
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Figure  A4.  Geometry  and  boundary  conditions  of  the 
terrain-vehicle  model 


e-  vsXl 


Figure  A5.  Slip  velocity  of  track  at  distance  X  from  the 
center  of  geometry 
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NORMAL  STRESS  DISTRIBUTION  WITHOUT  TRACK  TENSION 


NORMAL  STRESS  DISTRIBUTION  WITH  TRACK  TENSION 


Figure  A7.  Effect  of  track  tension  on  normal  stress  distribution 
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Figure  A9.  Track  speeds  and  velocities  of  an  arbitrary  point 

of  the  hull 
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APPENDIX  B:  DESCRIPTION  OF  THE  FIELD  DIRECT  SHEAR  DEVICE 


Background 

1.  The  terrain-vehicle  interaction  model  described  in  Appendix  A 
required  six  soil  parameters  as  input  that  have  to  be  determined 
experimentally .  The  parameters  are  (Figures  A1  and  A2,  Appendix  A): 

G  ,  initial  shear  stiffness  coefficient  (assumed  to  be  indepen¬ 
dent  of  rate  of  deformation) 

A  ,  material  constant  describing  the  maximum  shearing  strength  of 
the  material  at  very  high  normal  loading  (Equation  Al) 

M  ,  material  constant  related  to  the  parameter  A  and  to  the 
static  soil  cohesion  C  as  M  =  A  -  C 

N  ,  material  parameter  which  appears  in  Equation  Al 

C^  ,  increase  in  soil  cohesion  due  to  dynamic  loading  (maximum 
value  achieved  for  loading  rates  of  interest) 

A  ,  material  constant  describing  the  effects  of  rate  of  deforma¬ 
tion  on  the  cohesive  strength 

The  soil  parameters  G  ,  A  ,  M  ,  and  N  can  be  determined  from 

various  existing  laboratory  test  devices,  such  as  the  triaxial  shear 

device  or  direct  shear  device.  The  triaxial  shear  and  direct  shear 

devices,  however,  may  not  yield  the  same  values  of  G  ,  A  ,  M  ,  and 

N  for  identical  specimens  because  of  differences  in  test  boundary 

conditions.  The  stress  boundary  conditions  associated  with  the  direct 

shear  test  more  closely  approximate  the  stress  conditions  experienced 

by  the  soil  during  steering  of  track- laying  vehicles.  It  is,  therefore, 

more  appropriate  to  determine  these  parameters  from  direct  shear  tests. 

The  parameters  C,  and  A  can  only  be  determined  from  special  static 
d 

and  dynamic  triaxial  shear  tests  since  dynamic  direct  shear  devices  are 
not  presently  available.  Therefore,  to  adequately  determine  the  five 
soil  parameters,  two  separate  test  series  may  be  required: 
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a.  Direct  shear  tests  to  define  G  ,  A  ,  M  ,  and  N  . 

_b.  Static  and  dynamic  triaxial  shear  tests  to  define  Cj 
and  A  . 

It  should  be  noted  that  in  determining  and  A  from  triaxial  tests 

rather  than  direct  shear  tests,  it  is  assumed  that  these  parameters  are 
not  sensitive  to  test  boundary  conditions.  The  validity  of  this 
assumption  should,  of  course,  be  evaluated. 

2.  The  most  important  consideration  in  conducting  laboratory  soil 
tests  is  that  the  undisturbed  specimens  be  representative  of  the 
materials  over  which  the  vehicle  must  travel.  This  fact  implies  that 
the  upper  several  inches  of  surface  material  must  be  sampled,  trimmed 
to  necessary  specimen  size,  and  tested  in  the  laboratory.  Water  con¬ 
tent,  soil  structure,  density,  and  vegetation  root  systems,  all  of 
which  affect  material  response,  must  be  preserved.  With  this  in  mind, 
a  field-operated  direct  shear  device  capable  of  testing  a  variety  of  in 
situ  surface  soils  for  normal  loads  of  interest  was  designed  and 
fabricated.  The  description  of  the  device  and  the  procedure  by  which 
the  soil  parameters  can  be  determined  are  documented  in  this  appendix. 

Direct  Shear  Device 


Design  consideration 

3.  Previously  proposed  field  devices  were  considered  but  rejected 
because  of  one  or  more  of  the  following  reasons:  (a)  some  of  the  soil 
parameters  could  not  be  measured  and  hence  required  additional  tests; 

(b)  the  necessary  support  equipment  was  too  massive  to  be  easily  field 
transportable;  or  (c)  specimen  disturbances  were  encountered  before 
testing.  The  idea  of  creating  a  new  type  of  test  was  also  rejected 
because  any  new  device  would  contain  inherent  boundary  problems,  all  of 
which  would  have  to  be  evaluated  with  time  and  usage.  The  direct  shear 
device,  on  the  other  hand,  has  been  used  extensively,  and  it  is  a 
fairly  simple  test  to  run.  Furthermore,  the  four  basic  soil  parameters 
(G  ,  A  ,  M  ,  and  N)  could  be  measured  rather  directly  from  this 
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test.  Figure  R1  shows  a  sketch  of  the  field  device  that  was  fabricated 
as  a  result  cf  this  project.  Figure  B2  represents  photographs  taken  of 
the  device  during  the  conduct  of  actual  field  tests. 

Specimen  container 

4.  Specimen  configuration  was  the  first  consideration  made  in  the 
design  of  the  device.  It  was  assumed  that  in  many  cases  the  in  situ 
soil  could  not  be  sampled  without  disturbance;  therefore,  the  specimen 
container  would  have  to  be  placed  around  the  soil.  A  round  ring 
similar  to  a  coring  device  would  afford  the  least  chance  of  soil 
disturbance.  However,  the  stress  distribution  along  a  horizontal  plane 
of  a  circular  specimen  is  not  uniform.  To  reduce  the  nonuniformity,  a 
square-shaped  specimen  container  was  selected. 

5.  A  4-  by  4-in.  box  was  selected  in  order  to  keep  the  shear  and 
normal  loads  within  limits  of  interest  to  the  analysis  of  track-laying 
vehicles  and  at  the  same  time  retain  a  reasonably  large  specimen  size. 
The  use  of  deadweights  is  the  simplest  way  to  produce  normal  load,  but 
use  of  more  than  200  lb  in  weights  is  awkward  for  field  testing. 
Therefore,  with  the  weight  requirement  below  the  200-lb  limit,  normal 
stress  of  up  to  3.2  psi  can  be  produced  on  a  4-  by  4-in.,  r  16-sq- 
in.,  specimen.  However,  the  largest  particle  or  grain  size  permissible 
with  a  4-  by  4-in.  specimen  is  probably  1/2  in. ,  which  is  a  reasonable 
limit  for  most  terrains  of  interest. 

6.  The  overall  specimen  height  was  controlled  by  the  depth  of  the 
desired  shear  plane  as  directed  by  grouser  depth  ranging  from  approxi¬ 
mately  3/4  to  1-1/2  in.  The  compressibility  of  soil  could  significantly 
alter  this  depth,  but  for  estimation  purposes  the  depth  was  assumed  to 
be  no  greater  than  2  in.  Therefore,  the  height  of  the  upper  box  portion 
was  set  a*-  2  in.,  permitting  testing  of  depths  from  approximately 

1/4  to  2  in.  This  height,  of  course,  can  be  altered  should  particular 
site  conditions  dictate.  The  lower  box  portion  was  set  at  1-1/4-in. 
height,  including  the  cutting  edge.  A  1/8-in.  wall  thickness  was  used 
for  both  boxes. 

7.  vigure  83  presents  a  series  of  sketches  of  the  specimen  con¬ 
tainer  showing  the  various  stages  of  placement.  To  minimize  specimen 
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disturbance,  it  was  decided  to  use  the  specimen  cutting  box  as  the 
device  container  rather  than  remove  the  cutting  and  place  a  container 
over  the  specimen.  The  box  consists  of  three  parts:  (a)  a  lower 
portion  wi  .h  knife-sharp  edges  to  aid  in  cutting  the  soil,  (b)  an  upper 
portion,  and  (c)  an  outer  holder  to  keep  the  lower  and  upper  portions 
in  alignment.  The  box  is  alternately  pushed  and  trimmed  into  the  soil 
to  the  desired  depth.  Once  in  place,  the  outer  holder  can  vi  carefully 
removed,  leaving  the  two  boxes  on  the  specimen  with  the  joint  between 
the  lower  and  upper  box  portion  forming  the  shear  test  plane. 

Base 

8.  A  relatively  narrow  1-in. -thick  aluminum  plate  was  used  to 
construct  the  base  with  a  square  hole  at  one  end  to  fit  around  the  4- 
by  4-in.  lower  specimen  containers  (Figure  Bl) .  The  shear  loader  was 
attached  to  the  other  end  of  the  plat?.  A  second  1-in. -thick  aluminum 
yoke  was  constructed  to  fit  over  the  upper  specimen  container.  Set 
screws  through  the  yoke  serve  to  raise  the  yoke  off  the  baseplate,  thus 
minimizing  friction  between  the  surfaces.  The  shear  loader  attached  by 
cable  to  the  yoke  pulls  the  upper  specimen  while  the  base  reacts  against 
the  lower  specimen  container.  Guide  rails  along  the  edge  of  the  base 
ensure  that  no  torsional  shear  deformation  or  twisting  is  applied  to 
the  specimen. 

Shear  loader 

9.  An  electric  12-volt  boat  winch  was  incorporated  into  the  base 
as  the  shear  loader.  This  approach  is  the  simplest  for  providing  a 
shear  loader.  (If  necessary,  the  winch  can  be  replaced  with  a  more 
sophisticated  loader  custom-built  for  this  device.)  Currently,  the 
winch  is  capable  of  pulling  loads  up  to  2000  lb.  Static  loading  rates 
can  be  i pplied  by  manually  turning  the  winch  via  a  socket-rachet 
arrangement.  Fast  loading  rates  (approximately  300-600  msec  time  to 
peak  load  that  is  equivalent  to  a  strain  rate  of  0.5  to  1.5  per  sec, 
which  is  compatible  with  the  strain  rate  under  the  track)  can  be  applied 
using  the  electric  feature  of  the  winch.  The  power  is  supplied  ty  a 
12-volt  car  battery,  which  is  also  used  as  the  instrumentation  power 
supply . 
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Instrumentation 


10.  A  2-in.  travel  film  potentiometer  is  attached  to  the  base  and 
records  relative  movement  between  the  upper  specimen  holder  and  the 
base.  A  strain  gage  load  cell  attaching  the  winch  cable  to  the  specimen 
yoke  is  used  to  measure  shear  load.  A  compact,  two-channel  DC  instru¬ 
mentation  amplifier  is  used  for  signal  conditioning.  Output  is  recorded 
in  the  form  of  a  shear  load  versus  deflection  plot  on  a  commercially 
available  DC-operated  X-Y  plotter.  As  previously  mentioned,  a  simple 
car  battery  is  the  main  power  supply.  All  initial  testing  was  done  by 
recording  the  data  on  a  time  base  light  beam  strip  chart.  This  record¬ 
ing  procedure  was  later  dropped  since  the  loading  times  remained  fairly 
constant  on  the  soils  tested.  A  time  base  can  be  added  at  a  later  date 
through  the  use  of  a  frequency  oscillatory  and  an  X-Y-Z  recorder. 

Normal  load 

11.  A  series  of  steel  weights,  the  largest  weighing  57  lb  and  the 
smallest  weighing  8-1/2  lb,  was  fabricated  for  use  with  the  device. 

Guide  holes  and  studs  permit  stacking  and  centering  of  the  weights  on 
the  specimen  top  surface.  Although  a  variety  of  load  combinations  are 
possible,  most  tests  have  been  conducted  using  weights  totalling 
approximately  8.6,  36.6,  65.7,  and  122.7  lb  (i.e.,  normal  stress  levels 
of  0.54,  2.29,  4.11,  and  7.67  psi) . 

Measurement  of  Soil  Parameters 


12.  A  series  of  two  or  more  tests  is  required  at  a  site  to  define 
the  necessary  soil  parameters.  A  typical  testing  program  may  call  for 
the  conduct  of  four  fast  and  four  slow  tests  at  normal  stresses  of 
0.54,  2.29,  4.11,  and  7.67  psi.  For  each  test  an  X-Y  data  record  of 
shear  load  versus  deflection  is  obtained.  In  addition,  measurements  of 
soil  density  and  water  content  are  mace  on  each  test  specinen  (generally 
on  the  posttest  specimen  contained  in  the  upper  and  lower  specimen 
holders) . 
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13.  For  each  test  a  plot  of  shear  stress  versus  deflection  is 
obtained.  The  initial  slope  of  the  plot  defines  G  ,  the  peak  stress 
defines  the  maximum  shear  stress,  and  the  deflection  at  peak  stress 
divided  by  time  to  peak  stress  defines  the  deflection  rate.  A  table 
listing  of  each  test  is  used  to  summarize  the  data  and  contains  the 
specimen  number,  wet  density,  water  content,  dry  density,  normal 
load/stress,  maximum  shear  load/stress,  initial  G  ,  deflection  at  peak 
stress,  and  deflection  rate.  Figure  B4  presents  the  test  results 
obtained  from  the  series  of  field  tests  conducted  at  a  given  site. 

14.  The  analysis  plots  are  shown  graphically  in  Figure  B5.  A 

summary  plot  of  shear  stress  versus  shear  deformation  is  made  to 
obtain  either  static  or  dynamic  failure  envelopes.  From  the  static 
failure  envelope,  the  values  of  A  ,  M  ,  and  N  are  obtained.  The 
value  of  and  A  are  obtained  from  the  dynamic  failure  envelope  as 

shown  in  the  plot  of  cohesion  versus  rate  of  deformation.  The  value  of 
G  is  the  initial  slope  of  shear  stress-shear  deformation  curve  (Figure 
B4) . 


Conclusions 

15.  A  new  agility  model  for  track-laying  vehicles  was  developed 
that  required  soil  parameter  input  not  commonly  obtained  during  mobility 
studies.  It  v/as  realized  that  field  sampling  and  laboratory  testing 
would  not  always  be  possible  because  of  the  nature  of  very  near-surface 
soil  deposits.  An  approach  was  taken  to  use  a  conventionally  accepted 
test  to  define  the  parameters.  A  field-operable  direct  shear  device 
and  necessary  instrumentation  were  built  at  the  WES.  The  equipment  is 
fairly  compact  (can  easily  fit  into  a  car  or  truck),  operates  off  of  a 
car  battery,  is  relatively  simple  to  use,  saves  operation  time  compared 
with  comparable  laboratory  tests,  and  is  capable  of  directly  accessing 
the  soil  parameters. 

16.  The  device  has  been  used  to  conduct  some  123  tests  at  several 
different  locations.  Some  of  the  near-surface  soils  encountered  at 
these  locations  were  very  friable  and  impossible  to  sample  and  test  in 
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the  laboratory  by  conventional  means  without  excessive  disturbance. 
However,  the  direct  shear  device  performed  quite  well,  only  requiring 
extra  care  by  the  test  operator  not  to  disturb  the  specimen  during 
placement  of  the  device  base  over  the  sample  box  containing  the  soil 
specimen.  The  time  required  to  perform  a  test  was  approximately  one- 
half  hour.  This  time  is  especially  reasonable  when  compared  with  the 
time  it  takes  to  prepare  and  test  a  comparable  sample  in  the  laboratory, 
excluding  the  time  spent  obtaining  the  sample  in  the  field. 
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ORIGINAL 

GROUND 

SURFACE 


b.  Specimen  immediately  following  test  with  normal 
loads  removed 

Figure  112.  Photographs  of  the  direct  shear  device  taken 
during  actual  field  testing 
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OUTER  HOLDER 


CUTTING  EDGE 


Ni 


UPPER  SPECIMEN  CONTAINER 


LOWER  SPECIMEN  CONTAINER 


INTERFACE 
BETWEEN  THE 
UPPER  AND 
LOWER  SPECIMEN 
CONTAINERS 


mm 


a.  CONTAINER  WITH  OUTER  HOLDER  ON  SOIL  SURFACE 


b.  CONTAINER  DURING  PLACEMENT  -  ALTERNATELY  PUSHED  AND 
EXCESS  MATERIAL  SPOONED  (OR  TRIMMED)  OUT 


c.  CONTAINER  AT  DESIRED  DEPTH  -  OUTER  HOLDER  REMOVED  AND 
READY  FOR  PLACEMENT  OF  DIRECT  SHEAR  DEVICE 


Figure  B3.  Cross  sections  through  the  specimen  container 
showing  various  stages  of  piacement 
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DATA  FROM  A  SERIES  OF  FAST  TESTS  DATA  FROM  A  SERIES  OF  SLOW  TESTS 


TABLE  SUMMARIZING  PERTINENT  DATA  FROM  TESTS  CONDUCTED 
AT  A  GIVEN  SITE 


Figure  B4.  Presentation  of  test  results  obtained  from  the  direct  shear  device 


COHESION  MAXIMUM  SHEAR  STRESS,  SHEAR  STRESS 


LEGEND 
0  SLOW  TEST 
0  FAST  TEST 

REPRESENTATIVE  Y 
SELECTED  FROM 
DATA  TABLE 


RECOMMENDED  SOIL  PARAMETER 
FOR  A  GIVEN  SITE 

INITIAL  SHEAR 

STIFFNESS  G,  PSI/IN 

PARAMETER  A,  PS  I 

PARAMETER  M,  PSI 

PARAMETER  N,  1/PSI 

ADDED  DYNAMIC 
COHESION  Cd,  PSI 


Figure  B5.  Graphical  presentation  of  the  selection  of  recommended 
soil  parameters  based  on  the  field  data 
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APPENDIX  C:  SOIL  CLASSIFICATION  TEST  DATA 


1.  Appendix  C  contains  gradation  curves  and  soil  classification 
data  for  soil  samples  from  various  test  sites  (Figures  C1-C24) .  These 
test  sites  and  the  test  series  conducted  at  each  one  are  discussed  in 
Part  II  of  the  main  text. 


.  Gradation  curve  for  site  7A,  test  series  98-102, 


sample  1 
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Figure  C2. 


Gradation  curve  for  site  7A, 
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Figure  C3.  Gradation  curve  for  site  7A,  test  series  93-102,  sample  3 
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Figure  C4.  Gradation  curve  for  site  7B,  test  series  107-111,  sample  BA2 
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Figure  C5.  Gradation  curve  for  site  7B,  test  series  107-111,  sample  BA5 
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Figure  C7.  Gradation  curve  for  site  7A,  test  series  117-122,  sample  TP4 


Figure  C8.  Gradation  curve  for  site  7A,  test  series  117-122,  sample  TP2 
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Figure  C12.  Gradation  curve  for  site  7A,  test  series  131-134,  sample  STP2 
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Figure  C13.  Gradation  curve  for  site  7A,  test  series  131-134,  sample  STP4 


Figure  C14.  Gradation  curve  for  site  7B,  test  series  139-144,  sample  2A2 


Figure  C16.  Gradation  curve  for  site  7B,  test  series  139-144,  sample  2B3 
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Figure  C17.  Gradation  curve  for  site  7B,  test  series  139-144,  sample  1B1 
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Figure  C18.  Gradation  curve  for  site  7B,  test  series  13  iA4,  sample  CF1 
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Figure  C19.  Gradation  curve  for  site  7B,  test  series  139-144,  sample  CF4 
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Figure  C21. 
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Figure  C22 


Gradation  curve  for  site  8 


test  series  157-159,  sample  WES1 


Figure  C23.  Gradation  curve  for  site  8,  test  series  157-159,  sample  WES 5 
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Figure  C24 


Gradation  curve  for  site  8,  test  series  157-159,  sample  WES11 
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